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We  have  successfully  completed  the  study,  development 
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and  application  of  a  computer  controlled  Nc  r  pumped  LiNbOg 


tunable  source.  The  device  has  operated/i-eliably  for  more  than 
one  year.  It  tunes  over  a  1.4  -  O 


jrirr  range  and  operates  at 
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resolutions  to  less  than  0.08  -6m“  at  10  mj  output  at  10  pps. 


We  have  invented  and  demonstrated  a  new  optical  device, 

the  radial  birefringent  element.  We  have  invented  and  demonstrated 

a  computer  controlled  laser  attenuator. 

During  the  past  year  we  have  developed  the  theory  for  a 

slab  configuration  laser  source.  We  have  successfully  demonstrated 

the  slab  geometry  as  applied  to  Nd: Glass.  Work  is  continuing  on 
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the  understanding  and  implementation  of  this  high  energy,  high 
repet it ional  rate  Nd: Glass  slab  laser  source. 
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TUNABLE  OPTICAL  SOURCES 
Robert  L.  Byer 

FORWARD 

During  the  past  year  we  have  made  significant  progress 
in  three  aspects  of  tunable  laser  physics  :  operation  and 
measurements  with  a  NdrYAG  pumped  LiNbOg  Optical  Parametric 
Oscillator  Source  (0P0) ; invention,  development  and  demonstration 
of  a  new  optical  element  called  the  Radial  Birefringent  Element 
(RBE),  and  the  design  and  demonstration  of  a  slab  configuration 
Nd:Glass  laser  oscillator/amplifier  system. 

We  have  also  invented  and  developed  a  computer  controlled 
attenuator  and  have  transferred  that  technology  to  industry  through 
the  Stanford  University  License  and  Technology  Office.  We  have 
also  grown  hot  phasematching  LiNbOg  and  have  demonstrated  efficient 
second  harmonic  generation  of  a  cw  lamp  pumped  NdrYAG  laser  source. 

TUNABLE  OPTICAL  SOURCES 

A.  Statement  of  the  Problem  Studied 

The  work  during  the  past  year  aimed  at  using  nonlinear 
processes  to  generate  trunable  coherent  radiation.  The  work 
centered  on  the  pump  laser  sources  and  improvements  in* laser  source 
capability.  In  addition,  work  progressed  on  the  generation  of 
computer  controlled  coherent  tunable  radiation  via  the  NdrYAG 
pumped  LiNbOg  0P0. 


Work  began  on  the  study  of  Nd: Glass  as  a  tunable  source. 

The  slab  configuration  was  studied  with  a  design  goal  of  6  J 
per  pulse  at  10  pps  tunable  output  at  1.055.  ir  ±  100  cm”1. 

Stimulated  Raman  shifting  is  planned  for  frequency  extension  to 
the  infrared. 

LIST  OF  APPENDICES 

I.  M.  Endemann  and  R.L.  Byer,  "Remote  Atmospheric 

Measurements  in  the  Near  Infrared  Using  a  Reliable 
0P0  Source". 

IF.  Y.K.  Park,  R.L.  Byer,  R.S.  Feigelson  and  W.  Kway, 

"Growth  of  High  Phasematching  Temperature  LiNbOg 
Single  Crystals  and  Application  to  Efficient  Second 
Harmonic  Generation  of  NdrYAC". 

B.  Summary  of  Results 

During  the  past  year  we  completed  work  on  the  high  resolution 

operation  of  the  computer  controlled  Nd : YAG  pumped  LiNbOg  tunable 

source.  We  have  demonstrated  a  4500  cm-1  tuning  range  at  0.6  cm”1 

resolution;  a  0.08  cm-1  resolution  tunable  over  a  20  cm”1  range; 

more  than  one  year  of  operation  without  optical  damage  problems, 

% 

and  the  capability  to  make  precise  atmospheric  species,  temperature 
and  humidity  measurements.  The  results  are  summarized  in  Appendix  I 
and  have  been  presented  at  recent  conferences.  M.  Endemann' s  Ph.D. 
thesis,  in  preparation,  will  describe  the  tunable  source  in  detail. 
This  source  was  the  first  computer  controlled  tunable  device. 
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It  now  provides  the  widest  tuning  range  of  any  tunable  source 
and  has  operated  trouble-free  for  more  than  one  year.  A  number 


of  LiNbOg  0P0  sources  are  in  operation  around  the  world  for  remote 
measurement  and  laser  chemistry  research  applications.  This  work 
is  now  completed.  The  four  years  of  research  and  development 
have  led  to  a  major  advance  in  tunable  source  science  and  techno1 ogy. 

Research  into  methods  of  improving  laser  sources  continued 
during  the  past  year.  We  invented  a  new. optical  component,  the 
radi&l  birefringent  element,  and  applied  it  to  laser  resonator  design. 
We  invented  a  new  computer  controlled  laser  attenuator  and  demonstrated 
its  capability. 

We  studied  the  growth  of  hot  phasematching  LiNbOg  for  SHG 
of  Nd:YAG.  Our  work  led  to  the  solution  a  decade  old  growth 
problem  and  to  the  successful  SHG  of  a  Nd : YAG  laser  at  40%  efficiency. 
This  work  is  summarized  in  Appendix  II. 

Our  major  area  of  emphasis  was  the  theoretical  study  of  the 
slab  configuration  Nd  .’Glass  laser  source.  We  have  completed  a 
theoretical  description  of  the  slab  laser  and  have  developed  a 
first  generation  numerical  program  for  design  calculations.  We 
have  designed  and  constructed  a  slab  glass  oscillator.  The 
oscillator  verified  the  theoretical  design  calculations.  It 
clearly  showed  the  advantages  of  the  slab  geometry  for  glass  as 
well  as  for  high  average  power  Nd:YAG  and  other  crystalline  laser 
mat erials . 

We  are  actively  continuing  the  experimental  and  theoretical 
program  on  the  slab  geometry  laser  approach.  We  expect  to  complete 
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experimental  demonstration  of  the  slab  glass  system  in  the  near 
future  using  high  quality,  optimum  finished  LHG-5  Hoya  glass  slabs. 
This  work  will  then  be  written  as  a  theoretical/experimental  paper 
on  the  slab  geometry  laser  approach. 
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Remote  single-ended  measurements  of  atmospheric  temperature 
and  humidity  at  1.77  pm  using  a  continuously  tunable  source 


M.  Endemann  and  R.  L.  Byer 

Applied  Physics  Department,  Edward  L.  Ginzton  Laboratory,  Stanford  University,  Stanford.  California  94305 


Received  )une  9, 1980 

Simultaneous  remote  measurements  of  temperature  and  humidity  using  a  narrow-linewidth,  continuously  tunable, 
LiNbOa  optical  parametric  oscillator  as  a  transmitter  source  are  reported.  Relative  measurement  errors  of  1.0°C 
for  temperature  and  better  than  1%  for  humidity  over  a  45-sec  averaging  time  are  observed.  The  absolute  accuracy 
is  limited  by  the  accuracy  of  available  spectroscopic  data. 


This  Letter  describes  the  application  of  our  widely 
tunable  IR  air-pollution  monitoring  system1-2  to  the 
simultaneous  measurement  of  atmospheric  temperature 
and  humidity.  Simultaneous  remote  measurement  of 
these  two  important  meteorological  parameters-prom- 
ises  to  be  a  useful  tool  for  atmospheric  research. 

A  number  of  groups  have  proposed  and  analyzed 
methods  of  making  remote  temperature  measurements, 
and  preliminary  remote  temperature  measurements 
have  been  reported  using  Raman  lidar3  and  a  two-fre¬ 
quency  absorption  measurement  have  been  reported.4 
Remote  temperature  measurements  using  atmo¬ 
spheric-transmittance  measurements  of  three  wave¬ 
lengths  have  been  analyzed  by  Mason5  and  by 
Schwemmer  and  Wilkerson.6  Two  wavelengths  and 
i>2  are  selected  to  coincide  with  two  absorption  lines  of 
one  species  with  different  ground-state  energy  E°.  A 
third  wavelength  e0ff  is  tuned  away  from  the  absorption 
lines  and  is  used  to  normalize  the  other  two  measure¬ 
ments.  The  absorbance  A,  is  given  by 

Ai(T)  =  NL  odT)  =  ln(S0ff)  -  ln(S,),  (1) 

where 

S,  is  the  signal  received  at  wavelength  t\, 

N  is  the  density  of  the  species, 

L  is  the  pathlength,  and 
«r,(T)  is  the  absorption  cross  section  at  and 
temperature  T. 

The  temperature  dependence  of  the  absorption  cross 
section  is  given  by  the  Boltzmann  distribution  and  a 
partition  function  Q(T)  as 

<r,(T)  =  <n(T0)Q(T)expl(E\/k)(l/T0  -  1/T ))  (2) 

for  an  arbitrary  reference  temperature  To,  where  E '  is 
the  lower-state  energy  of  the  transition  at  e,  and  k  is  the 
Boltzmann  constant. 

From  the  two  measured  values  A,  (T),  we  can  deter¬ 
mine  the  temperature  T  as 


To 


,  feTp  ./AilT)  ■  ffi(To) 

E]  -  E\  lnU2(T)j  lno2(T0). 


(3) 


or,  for  (T  -  To)/T0  «  1, 


T  =  CT  +  D,  (4) 

with  7=  In  Ai(T)  -  In  A2(T),  C  =  kTl/(E\  -  E\),  and 
D  =  C[ln  ff2(T0)  —  In  <Ti( T'o) J  +  To(fe). 

As  expected,  the  accuracy  of  the  temperature  mea¬ 
surement  is  inversely  proportional  to  the  difference  in 
lower-state  energies,  E\  -  E%.  Once  the  temperature 
is  known,  it  is  straightforward  to  calculate  the  density 
of  the  species. 

Using  selected  H20  transitions  permits  the  simul¬ 
taneous  measurement  of  temperature  and  H20  con¬ 
centration  and  thus  of  relative  humidity.  Because  of 
the  complexity  of  the  spectrum,  we  used  tapes  supplied 
by  the  U.S.  Air  Force  Cambridge  Research  Laboratory 
(AFCRL)  to  select  H20  lines  with  appropriate  absorp¬ 
tion  cross-section  and  ground-state  energy  differences. 
For  our  measurements,  the  optimum  absorption  lines 
were  found  in  the  1.9 -/rm-band  region.  The  J  =  11  line 
at  5651.33  cm-1,  with  a  lower-state  energy  of  1999  cm-1, 
and  a  nearby  <7  =  6  line  at  5650.41  cm-1,  with  a  543  cm-1 
lower -state  energy,  were  selected  since  their  absorption 
cross  sections  at  room  temperature  are  of  similar 
strength  and  their  absorbance  is  approximately  unity 
at  average  H20  concentrations  over  a  1  -km  path.  The 
1456-cm_l  difference  in  lower-state  energy  yields  a 
0.8°C  temperature  error  for  a  1%  accuracy  in  deter¬ 
mining  the  absorbance  at  and  i*j>. 

One  problem  does  arise  with  thi?  selection  of  lines. 
At  each  wavelength,  not  only  do  the  setecved  lines  con¬ 
tribute  to  the  measured  absorbance;  the  wings  from 
some  strong  nearby  lines  and  many  weaker,  but  not 
negligible,  lines  also  contribute.  We  used  the  spectro¬ 
scopic  data  contained  in  AFCRL  tapes7  in  a  computer 
program  to  simulate  this  interference.  Figure  1  shows 
a  simulation  of  a  temperature  measurement  using  the 
same  three  wavelengths  as  were  used  in  the  actual  ex¬ 
periment.  Here  T  (see  Eq.  (4)]  is  plotted  against  the 
actual  temperature  T.  Overlapping  of  the  two  main 
absorption  lines  and  the  influence  of  other  lines  and 
second-order  effects  make  7’  a  nonlinear  function  of  T. 
However,  the  linear  approximation  Eq.  (4)  holds  to 
better  than  1  °  accuracy  from  -20  to  30°C.  Variat  ion 
of  the  spect  ral  resolut  ion  in  the  computer  simulation 
shows  only  small  changes  of  the  constants  C  and  D  for 
temperat  ures  almve  — 20°C.  Spectroscopic  data  used 
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Fig.  1.  Temperature  calibration  using  spectroscopic  data. 
For  only  two  absorption  lines,  the  temperature  dependence 
is  T  =  CT  +  D.  The  interference  of  other  lines  and  the  sec¬ 
ond-order  effects  of  the  approximation  lead  to  the  slightly 
curved  calibration  line.  Different  probing  resolutions  change 
the  calibration,  but  temperature  shifts  are  small  above 
— 20®C. 


in  this  simulation  have  line-strength  errors  of  about 
20%,  whereas  the  lower-state  energies  are  known  to  0.02 
cm-1.  The  line-strength  error  affects  the  constant  D, 
whereas  the  value  for  the  proportionality  constant  C  is 
accurate.  In  actual  measurement,  we  had  to  use  a  value 
for  D  that  is  20°C  less  than  that  predicted  in  the  simu¬ 
lation.  We  are  working  to  improve  the  spectroscopic 
data  by  using  our  tunable  source. 

The  experimental  setup  is  an  improved  version  of  the 
LiNbOa  optical  parametric-oscillator  (OPO)-based 
air-pollution  monitoring  system.1-2  An  unstable  reso- 
nator-Nd:YAG  laser  is  used  to  pump  an  angle  phase- 
matched  LiNbO-t  OPO.  The  OPO  is  continuously 
tunable  over  a  4500-cm-l-wide  range,  from  1.42  to  4.0 
Mm.  It  generates  10-mJ  output  energy  in  a  15-nsec-long 
pulse  with  a  10-Hz  repetition  rate.  It  has  operated  re¬ 
liably  for  more  than  one  year.  Details  of  the  tunable 
source  have  been  reported  earlier.2-8 

The  linewidth  of  the  OPO  is  narrowed  to  0.7  cm-1  by 
using  an  intracavity  grating.  A  tilted  etalon  further 
reduces  the  output  linewidth  to  0.1  cm-1.  This  source 
is  continuously  tunable  by  using  a  PDP11  E10  mini¬ 
computer  to  set  crystal  angle,  grating,  and  6 talon.  The 
OPO  output  is  expanded  to  a  2.5-cm-diameter  beam 
and  directed  to  the  roof  of  our  laboratory,  where  it  is 
coaxially  transmitted  through  the  telescope.  The 
40-cm-diameter  receiving  telescope  focuses  the  received 
light  onto  a  liquid-nitrogen-cooled  InSb  detector. 

Since  the  output  energy  of  the  OPO  is  not  adequate 
for  differential  absorption  measurements  using  Mie 
backscattering,  we  were  limited  to  long-path  absorption 
measurements  using  a  building  at  a  distance  of  775  m 
as  a  backscattering  target.  The  height  above  ground 
level  of  this  path  varies  from  10  to  60  m, 

A  computer  was  used  to  tune  the  OPO  and  for  data 
taking,  averaging,  storage,  and  display.  An  extensive 
interactive  routine  was  developed  to  obtain  atmo¬ 
spheric-transmission  spectra  and  to  make  automatic 
temperature  and  humidity  measurements.  The  lower 


trace  in  Fig.  2  shows  an  example  of  an  atmospheric- 
transmission  spectrum  with  a  0.1  -cm  resolution  tuned 
over  a  17-cm-1  range.  The  upper  trace  is  a  simulated 
transmission  spectrum  made  by  using  the  AFCRL 
tapes.7  We  used  scans  similar  to  this  to  calibrate  the 
tuning  of  the  OPO  and  to  identify  the  temperature- 
sensitive  lines. 

Once  the  wavelengths  of  these  lines  are  determined 
and  the  OPO  aligned,  temperature  and  humidity  mea¬ 
surements  are  made  under  computer  control.  The  re¬ 
turn  signal  at  each  frequency  is  averaged  over  several 
shots  before  we  tune  to  another  line.  The  time  of  this 
average  has  been  optimized  to  minimize  errors  that  are 
due  to  water-vapor  fluctuations  but  to  achieve  maxi¬ 
mum  time  resolution  of  the  measurement.  The  OPO 
is  now  tuned  after  30  shots  or  3  sec,  which  is  a  good 
compromise  in  calm  atmospheric  conditions.  After  all 
three  lines  are  probed,  temperature  and  humidity  are 
calculated,  displayed,  and  stored.  Further  averaging 
can  be  done  later;  this  reduces  the  storage  space  re¬ 
quired.  Also,  the  standard  deviation  of  each  quantity 
resulting  from  return-signal  fluctuations  is  determined 
and  stored. 

Figure  3  shows  the  record  of  temperature  and  hu¬ 
midity  from  the  morning  of  March  1 1 , 1980.  Each  point 
represents  150  shots /line  or  averaging  over  45  sec.  The 
standard  deviation  was  1.1°C  in  the  beginning  but  in¬ 
creased  to  1.4°C  at  0800  h.  The  outside  weather  was 
cloudy  at  first,  with  rain  starting  at  0740.  The  rain  then 
increased  until  it  was  so  strong  at  0810  that  the  mea¬ 
surement  had  to  be  interrupted.  At  0820  the  rain 
stopped,  the  cloud  cover  slowly  dissipated,  and  the  wind 
increased.  The  remotely  measured  temperature  is 
compared  with  a  record  of  a  thermograph  at  the  tele¬ 
scope.  The  general  agreement  is  good.  Some  dis¬ 
crepancies  in  the  two  measurements  are  expected  since 
the  remote  measurement  averages  over  a  path  from  10 
to  60  m  above  the  ground  and  775  m  in  length,  whereas 
the  thermograph  is  located  on  the  roof  of  a  building. 

The  calibration  of  our  data  is  done  by  using  a  large 
number  of  remotely  measured  temperature  values  (we 
have  accumulated  a  total  of  10  h  of  data)  and  comparing 


ATMOSPHERIC  TRANSMISSION  OVER  I.SSfcm  LONG  PATH 


S650  5655  5666 

FREQUENCY  UsC't 

Fig.  2.  High-resolution  atmospheric  scan  over  1.55-km  path. 
The  OPO  linewidth  is  0. 1  cm" 1 .  The  OPO  transmission  scan 
(lower  truce)  is  compared  with  a  simulation  using  AFCRL  file 
data  (upper  trace). 
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Fig.  3.  Temperature  and  absolute  humidity  measurement 
from  the  morning  of  March  1 1 , 1980.  Temperature  data  are 
compared  with  records  from  thermograph  at  receiving  tele¬ 
scope. 


them  with  the  thermograph  data.  We  found  that  the 
best  value  for  the  constant  D  is  35°C,  which  is  19°C  less 
than  the  value  predicted  from  the  available  spectro¬ 
scopic  data. 

Gbr  data  for  humidity  measurement  have  not  been 
checked  with  an  independent  instrument  so  far.  We 
used  the  absorption  strength  for  the  5650.41-cm"1  line 
from  the  AFCRL  tape,7  so  the  accuracy  is  estimated  to 
be  ±20%.  We  compared  our  values  with  relative  hu¬ 
midity  and  temperature  measurements  at  nearby  me¬ 
teorological  observation  stations.  Although  the  trend 
of  station  data  confirmed  ours,  the  local  fluctuations 
were  too  big  to  improve  the  calibration  for  our  data. 

We  have  demonstrated  the  feasibility  of  measuring 
temperature  and  humidity  simultaneously  by  using  a 
high-resolution,  continuously  tunable  source.  Ab¬ 
sorption  lines  in  the  1.9-pm  H20  band  have  large 
enough  energy  differences  in  the  lower  states  to  permit 
temperature  measurements  to  be  made  with  1°C  ac¬ 
curacy.  Fluctuations  of  the  H20  concentration  do  not 
cause  measurement  problems  as  long  as  the  time  for 


each  measurement  is  shorter  than  the  period  of  these 
fluctuations. 

Our  Nd:YAG-pumped,  LiNbOn  OPO-based,  remote 
monitoring  system  proved  to  be  reliable.  Operation  at 
narrow  linewidth  was  necessary  to  resolve  the  H2O  lines, 
and  the  use  of  a  continuously  tunable  source  was  es¬ 
sential  for  optimum  line  selection.  Finally,  full  com¬ 
puter  control  was  utilized  for  automatic,  long-term  re¬ 
mote  measurements. 

In  conclusion,  we  have  measured  temperature  to  1° 
accuracy  with  a  time  resolution  of  45  sec  and  deter¬ 
mined  the  H20  concentration  with  a  relative  error  of  less 
than  1%.  However,  the  absolute  accuracy  of  the  hu¬ 
midity  data  is  limited  at  present  by  the  accuracy  of  the 
available  spectroscopic  data. 

We  wish  to  acknowledge  the  support  of  the  U.S.  Army 
Research  Office  through  contract  #DAAG29-77-G- 
0181  and  the  Electric  Power  Research  Institute  for  the 
continued  use  of  equipment. 
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Computer-controllable  wedged-plate  optical  variable 
attenuator 


K.  Bennett  and  R.  L.  Byer 


A  computer-controllable  variable  optical  attenuator  has  been  developed  and  tested  that  has  a  tvide  dynamic 
range,  wide  spectral  range,  and  is  suitable  for  applications  with  high  peak  and  average  power  laser  sources. 
The  device  is  based  on  Fresnel  transmission  through  two  pairs  of  wedged  plates.  A  35-dB  dynamic  range, 
an  insertion  loss  of  1%,  a  precision  of  better  than  1%  and  beam  offset  and  deflection  of  <0.5  mm  and  0.5 
mrad,  respectively,  are  demonstrated. 


I.  Introduction 

The  ideal  optical  variable  attenuator  should  provide 
wide  dynamic  range,  wide  spectral  range,  resettability, 
and  high  accuracy,  yet  exhibit  a  simple  optical  and 
mechanical  design.  We  added  the  additional  require¬ 
ments  that  the  optical  variable  attenuator  operate  with 
high  peak  and  average  power  inputs  produced  by  Q- 
switched  laser  sources  and  not  significantly  alter  the 
beam  profile,  position,  or  direction. 

With  these  constraints  in  mind,  we  reviewed  possible 
methods  of  applying  the  three  fundamental  interactions 
that  can  be  used  to  attenuate  optical  radiation:  ab¬ 
sorption;  Fresnel  reflection  at  a  surface;  and  interfer¬ 
ence.  Attenuation  by  absorption  and  interference  were 
eliminated  as  not  being  consistent  with  the  require¬ 
ments  for  high  peak  and  average  power  and  wide  spec¬ 
tral  range. 

Fresnel  reflection  is  the  basis  for  many  forms  of  at¬ 
tenuators  including  those  based  on  birefringent  crystals 
and  polarization-dependent  total  internal  reflection.1 
We  chose  to  avoid  the  use  of  crystals  because  of  limited 
spectral  range  and  optical  quality.  Instead  we  con¬ 
centrated  on  Fresnel  reflection  near  the  critical  angle 
in  a  homogeneous  dielectric  medium. 

The  variable  gap  spacing  total  internal  reflection 
attenuator2  is  spectrally  sensitive  and  requires  tight 
fabrication  tolerances.  On  the  other  hand,  a  variable 
incident  angle  wedged-plate  attenuator  utilizing  Fresnel 
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reflection  meets  all  of  the  above  design  criteria  when 
properly  executed;  beam  deflection  and  offset  are  lim¬ 
ited  only  by  fabrication  errors. 

The  theory,  design,  and  performance  of  the 
wedged-plate  optical  attenuator  are  described  in  Secs. 
II-IV. 

II.  Theory 

The  power  transmittance  (Ts)  for  a  single  surface 
depends  upon  the  polarization,  the  ratio  of  the  indices 
of  refraction  n,  and  the  angle  of  incidence  0,.  The 
transmittance  for  light  polarized  parallel  and  perpen¬ 
dicular  to  the  plane  of  incidence  is3 


Tsi(6i,n)  »  — , 
sin2(0,  +  6, ) 

(1) 

T,  .fftnl—  Ts±(t>i<n) 

Ts,<<u)  co* He,-et)’ 

(2) 

where  Snell’s  law  in  the  form 

n  sin®,  =  sin®, 

(3) 

applies,  and  n  =  nt/nt  is  the  index  ratio  with  n,  and  n, 
the  refractive  indices  of  the  incident  and  transmitting 
media. 

Snell’s  law  leads  to  total  internal  reflection  when  no 
real  solution  for  ft,  exists  and  occurs  only  if  n,  <  n,  or 
when  n  <  1.  To  achieve  total  internal  reflection  we 
designed  the  attenuator  to  use  wedged  plates  of  wedge 
angle  <v. 

The  transmittance  of  a  single  wedged  plate  T„  is  the 
product  of  the  transmittance  of  two  surfaces, 

-  T.(#„n>Ts(»,  +  <t,  1/n),  (4) 

and  the  transmittance  of  the  entire  attenuator  T  is  the 
product  of  the  transmitlances  of  each  of  the  four  iden¬ 
tical  wedges: 
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Pig.  1.  Schematic  of  wedged-plate  attenuator  optical  layout.  Axes 
of  rotation  are  shown  near  base  of  incident  and  exit  prisms. 


T»  |T„(fli,n,  a)]4.  (5) 

We  define  the  critical  angle  of  incidence  df'1  for  the 
wedged  plate  as  the  minimum  value  of  0;,  which  results 
in  total  internal  reflection 

sinflr1 »  n  sin(u  -  a),  where  n  sinu  s  1.  (6) 

In  the  next  section  we  apply  the  above  results  to  the 
design  of  a  wedged-plate  attenuator. 


III.  Design 


A.  General  Design  Procedure 

Figure  1  shows  the  schematic  of  the  wedged-plate 
attenuator  design.  The  wedge  angle  «,  wedge  length 
L,  wedge  spacing  S,  wedge  overlap  D,  and  aperture  a 
required  to  accept  the  beam  are  shown.  The  aperture 
may  be  removed  from  the  problem  by  defining  nor¬ 
malized  coordinates  /  =  L/a,  s  =  S/a,  and  d  &  D/a.  The 
wedge  length  L’  must  be  sufficiently  longer  than  L  to 
account  for  wedge  thickness  and  fabrication.  To  design 
the  attenuator  the  four  parameters  a,  l,  s,  and  d  must 
be  chosen.  The  primary  considerations  are  the  maxi¬ 
mum  attenuation  required  and  the  allowable  interfer¬ 
ence.  For  large  aperture  rat  ios  (/  £  15),  nonuniformity 
of  the  attenuation  arising  from  beam  divergence  may 
be  a  consideration  at  extremely  high  attenuations. 
Insertion  loss  for  the  polarization  in  the  plane  of  inci¬ 
dence  is  <1%  for  l  >  4  and  n  *  1.5  and  is  relatively  in¬ 
dependent  of  the  wedge  angle  for  the  orthogonal  po¬ 
larization. 

The  fractional  beam  expansion  in  the  plane  of  inci¬ 
dence  is  given  by 


da 


a 


COS2flj 

cos30i, 


A»i, 


«in9i,  ■  n  »in(n  +  i>), 


(7) 

(8a) 


it  aitw  *  ainOj, 


(8b) 


where  0t  2  is  the  angle  between  the  normal  to  the  wedge 
and  the  beam  in  the  gap  between  wedges  in  a  pair,  and 
A 6i  is  the  full  angle  of  the  incident  laser  beam  diver¬ 
gence  (radians).  The  beam  expansion  A a/a  is  <2OOA0, 
for  the  designs  considered  here  and  is  not  an  important 
consideration. 

The  wedge  angle  a  is  chosen  so  that  the  laser  beam 
is  entirely  intercepted  at  the  critical  angle  given  by  Eq. 
(6).  The  wedge  angle  is 

a  •*  u  -  w,  (9») 

where 

sinw  “  1/n  »in(arcsee  /).  (9b) 

The  wedge  spacing  s  is  determined  from  the  maxi¬ 
mum  allowable  interference.  Interference  occurs  when 
beams  multiply  reflected  by  the  parallel  surfaces  on 
either  side  of  the  gap  are  not  offset  by  at  least  the  ap¬ 
erture  a.  Assuming  a  uniform  intensity  across  the  ap¬ 
erture  and  zero  elsewhere,  the  fraction  /  of  the  trans¬ 
mission  through  the  aperture,  which  is  doubly  reflected, 
is 

I m  r*(0»,)(l  -  2*  tanfft'Corfj)  (10) 

if  2s  tan0<2  cos0,  <  1  and  is  zero  otherwise.  Here  r  =  1 
-  T,(0,„n). 

After  choosing  a  maximum  value  of  /  for  each  polar¬ 
ization,  the  required  wedge  spacing  is  found  by  nu¬ 
merically  solving  for  the  maximum  value  of  s  as  a 
function  of  (>t2.  For  the  polarization  in  the  plane  of 
incidence,  s  has  a  single  maximum  for  0,2  greater  than 
Brewster’s  angle. 

A  minimum  wedge  overlap  d  is  required  in  each 
wedged-plate  pair  to  avoid  beam  clipping  by  the  second 
wedge.  Snell’s  law  and  geometry  yield 

d  m  l/confl,  -  *  tan 0«,  (11) 
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for  the  beam  walk  where  the  thickness  of  the  optics  has 
been  neglected.  This  function  has  a  single  maximum  ■ 
in  the  domain  of  interest  and  is  again  computed  nu¬ 
merically  using  the  value  of  s  determined  above. 

The  maximum  attenuation  is  limited  by  clipping  of 
the  outside  edge  of  the  second  wedge.  This  occurs 
for 

Unfl„  a  ((/  -  <0/4  (12) 

The  maximum  attenuation  is  found  by  applying  Eq.  (5) 
with  6i  calculated  from  Eqs.  (8a)  and  (8b). 


B.  Numerical  Design 

Figures  2  and  3  can  be  used  to  design  attenuator  op¬ 
tics.  These  were  computed  assuming  n  =  1.45  and  a  1% 
maximum  interference  in  the  Brewster  polarization, 
which  implies  a  3-6%  maximum  interference  in  the  or¬ 
thogonal  polarization.  The  procedure  is  as  follows: 

(1)  Determine  the  desired  maximum  attenuation 
and  find  the  required  normalized  length  from  Fig.  2, 
which  was  calculated  according  to  the  method  of  the 
previous  section,  specifically  Eqs.  (5)  and  (9)-(12).  The 
wedge  length  must  be  la  +  wedge  thickness  +  fabrica¬ 
tion  tolerance. 

(2)  Determine  the  wedge  angle  «  from  Eq.  (9). 
This  is  shown  along  the  top  axis  of  Fig.  3. 

(3)  Read  off  the  spacing  and  overlap  requirements 
from  Fig.  3. 

Interference  can  be  reduced  by  increasing  the  spacing. 
However,  an  increase  in  wedge  spacing  also  reduces  the 
maximum  attenuation. 

Using  these  equations  and  design  procedure,  we 
evaluated  the  attenuator  performance  as  a  function  of 
wedge  angle  and  index  of  refraction  for  a  number  of 
factors  of  interest  including  attenuator  dynamic  range, 
insertion  loss,  beam  displacement  and  angular  devia¬ 
tion,  and  acceptance  angle.  One  design  that  provided 
a  good  compromise  between  attenuation  range  and  ease 
of  optical  and  mechanical  fabrication  is  described 
next. 


N0RMALI2E0  WEDGE  LENGTH  (I) 

Fig.  2.  Maximum  attenuation  vs  aperture  ratio  /  for  n  =  1.45, 1.0% 
maximum  interference  in  Brewster  polarization  (see  text). 


WEDGE  ANGLE  o  (DEGREES) 


Fig.  3.  Normalized  wedge  spacing  s,  normalized  wedge  overlap  d, 
and  wedge  angle  «  for  the  specifications  of  Fig.  2. 


C.  Specific  Wedged-Plate  Attenuator  Design 
Figure  1  shows  a  schematic  of  the  wedged -plate  at¬ 
tenuator.  Four  identical  uncoated  fused  silica  wedges 
were  fabricated  with  A/10  quality  surfaces  and  wedge 
angles  «  =  1°.  This  design  is  not  strictly  according  to 
the  criteria  of  the  above  section,  because  the  attenuator 
was  designed  to  work  well  both  for  a  1-mm  diam,  1-mrad 
divergence  beam  and  6-mm  diam,  0.5-mrad  divergence 
beams  with  a  10-mm  positioning  aperture.  The  pairs 
of  wedges  were  mounted  on  counterrotating  shafts  lo¬ 
cated  to  minimize  beam  displacement  and  optical 
alignment  problems.  The  shafts  were  driven  by  inde¬ 
pendent  stepper  motors  geared  through  antibacklash 
worm  gears  to  a  ratio  of  2400  steps/rad. 


\ — APERTURE  - 
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TO  COMPUTER  AND 
MANUAL  CONTROLS 

.TO  AMPLIFIER 
AND  COMPUTER 

Fig.  4.  Schematic  of  optical  and  electrical  system  used  for  attenuator 
calibration  measurements. 
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Fig.  5.  Measured  transmittance  7V  =  0.6  sinWir/200)  vs  stepper  motor  position  N: 
•  ,  data  points;  — ,  theoretical  transmittance;  rms  error  is  0.008. 


The  wedges  were  fabricated  75  mm  in  length. 
Alignment  of  the  wedges  within  a  pair  must  be  within 
the  diffraction  angle  of  the  incident  beam  to  avoid  sig¬ 
nificant  beam  deflection.  The  wedges  do  induce  some 
beam  expansion  in  the  plane  of  incidence.  However, 
for  incident  light  polarized  perpendicular  to  the  plane 
of  incidence,  beam  expansion  is  <0.5  mm  for  a  1-mrad 
divergence  input  beam  at  the  38-dB  attenuation  level 
and  decreases  rapidly  at  lower  attenuations.  The  cal¬ 
culated  insertion  loss  is  24%  for  a  dynamic  range  of  35 
dB  in  this  polarization. 

For  a  beam  polarized  parallel  to  the  plane  of  incidence 
and  incident  at  Brewster’s  angle,  the  calculated  inser¬ 
tion  loss  of  the  attenuator  is  <0.1%.  The  attenuator 
provides  15-dB  dynamic  range  in  this  polarization. 


IV.  Measurement  Results 

The  performance  of  the  attenuator  was  verified  using 
a  polarized  He-Ne  laser  source  and  a  wide  dynamic 
range  silicon  diode  detector.  The  attenuator  was  con¬ 
trolled  hy  a  PDP11/10  minicomputer  that  set  the 
stepper  motors  and  recorded  the  transmitted  and  ref¬ 
erence  powers  at  1-sec  intervals  for  each  stepper  motor 
position.  A  schematic  of  the  experimental  setup  is 
shown  in  Fig.  4. 

The  attenuator  was  placed  in  the  He-Ne  beam  path 
in  approximate  alignment.  Final  angular  alignment  of 
the  prisms  was  accomplished  to  within  0.5  mrad  hy 
manually  stepping  the  motors  to  critical  angle  inde¬ 
pendently.  This  procedure  provided  an  absolute  or¬ 
ientation  for  the  prisms  and  proved  to  be  very  repeat- 
able. 

The  laser  output  was  then  blocked,  and  the  dark 
current  was  measured  for  the  signal  and  reference 
channels.  These  levels  were  used  ns  the  zero  signal 
levels  for  the  subsequent  measurements. 


To  show  the  versatility  and  the  dynamic  range  of  the 
attenuator,  two  driving  patterns  were  chosen.  The  first 
in  a  sine  attenuation  curve,  and  the  second  is  a  scan  at 
constant  angular  rate.  The  sine  attenuation  routine 
calculates  the  proper  number  of  motor  steps  to  generate 
a  power  transmittance  function  Tjv  =  0.6  sin(N7r/200) 
accurate  to  ±0.0025.  The  data  obtained  by  this  routine 
are  shown  in  Fig.  5  for  polarization  perpendicular  to  the 
plane  of  incidence.  The  theoretical  curve,  which  is 
drawn  as  a  solid  line  in  Fig.  5,  is  calculated  according  to 
the  actual  motor  setting  and  scaled  to  account  for  a 
measured  1%  scatter  loss.  The  average  rms  error  be¬ 
tween  the  actual  and  calculated  attenuation  is 

_  rih«n,)2j1/2 

V  n-  f  1  -0-008' 

The  error  agrees  well  with  the  average  standard  devia¬ 
tion  of  the  data  of  0.01  and  the  expected  average  digi¬ 
tization  uncertainty  of  0.005. 

Figure  6  shows  the  attenuation  obtained  for  polar¬ 
ization  perpendicular  to  the  plane  of  incidence  as  a 
function  of  scan  angle  relative  to  the  critical  angle.  For 
the  fused  silica  wedges  used,  n  =  1.45  and  0""  =  79.0° 
for  a  =  1°.  The  attenuation  is  plotted  on  a  semilog 
scale  to  illustrate  the  wide  dynamic  range.  To  retain 
signal  levels  for  accurate  digitization,  the  amplifier  gain 
settings  were  changed  and  normalization  chosen  for 
each  setting  to  match  the  data  points  in  overlapping 
regions.  The  maximum  attenuation  was  ~4000:1,  and 
the  insertion  loss  for  this  polarizat  ion  was  25%,  which 
yields  a  dynamic  range  of  35  dB.  The  average  rms  error 
( was  7, 1.5, 2,  and  0.7%  for  the  >30-,  20-30-,  10-20-,  and 
10-dB  attenuation  ranges.  The  error  is  larger  in  the 
10-20-dB  range,  because  the  amplified  signal  was 
smaller,  which  led  to  increased  electrical  noise  and  di¬ 
gitization  errors.  Data  in  the  15-20-dB  range  were 
taken  at  a  higher  amplifier  gain  setting  and  had  an  rms 
error  of  1%.  All  averages  are  within  one  standard  de¬ 
viation  of  the  theoretical  values. 
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Fie.  6.  Attenuation  in  dB  vs  angle  from  critical  angle  for  incident  He-  Ne  beam  polarized  per¬ 
pendicular  to  the  plane  of  incidence:  •,  measured  data  points;  — ,  theoretical  attenuation. 


In  addition,  all  points  in  the  >30-dB  range  showed 
a  transmittance  in  excess  of  the  theoretical  value.  This 
is  due  to  some  residual  component  of  the  orthogonal 
polarization.  For  example,  an  orthogonal  polarization 
component  of  only  one  part  in  seven  hundred  would 
account  for  the  7%  extra  transmittance  in  the  >30-dB 
range. 

Similar  data  were  taken  for  the  Brewster  polarization, 
and  a  maximum  attenuation  of  15  dB  was  measured  as 
expected.  The  measured  insertion  loss  was  near  1%, 
which  is  significantly  greater  than  the  expected  0.1% 
loss.  However,  eight  wedge  surfaces  provide  ample 
opportunity  for  scatter  losses  to  contribute  to  the  at¬ 
tenuation. 

The  attenuator  performance  was  also  verified  with 
a  Q-switched  Nd:YAG  laser  source  operating  at  70 
mJ/pulse  in  a  6-mm  diam  beam.  Since  only  uncoated 
fused  silica  surfaces  are  employed,  the  attenuator 
should  operate  at  incident  energy  densities  limited  only 
by  the  surface  damage  threshold  of  fused  silica  >3  J/ 
cm2. 

V.  Conclusion 

We  developed  a  simple  but  versatile  computer-con¬ 
trollable  optical  variable  attenuator  suitable  for  use  with 
virtually  all  collimated  light  sources,  including  high 
power  Q-switched  lasers.  The  spectral  response  is  ex¬ 
tremely  flat  and  can  he  accurately  calculated  from  the 
refractive  index.  The  useful  spectral  range  is  very 
broad  limited  only  by  the  transmittance  range  of  the 
wedge  material.  The  optical  fabrication  and  mechan¬ 
ical  design  are  straight  forward,  and  the  attenuator  docs 
not  require  dielectric  or  metallic  coatings.  Fresnel’s 
equations  and  a  knowledge  of  the  wedged-plate  indices 
of  refraction  lead  to  an  absolute  value  of  attenuation 
with  an  accuracy  limited  by  the  knowledge  of  incident 
polarization.  Transmittances  precise  to  <2%  are  easily 
obtained.  Absolute  accuracy  can  l>e  improved  by  using 
two  attenuators  in  tandem  with  relative  planes  of  inci¬ 
dence  rotated  by  precisely  90°  to  eliminate  polarization 
dependence. 


The  attenuation  has  a  scatter-limited  insertion  loss 
for  the  Brewster  polarization  of  <1%  and  a  dynamic 
range  of  15  dB.  The  non-Brewster  polarization  inser¬ 
tion  loss  is  25%  with  a  dynamic  range  of  35  dB.  Beam 
deflection  and  offset  are  limited  by  fabrication  and 
optical  alignment  tolerances.  A  simple  hand  assembly 
technique  resulted  in  negligible  beam  offset  and  de¬ 
flection  except  at  the  highest  attenuation  levels  where 
the  measured  offset  is  <0.5  mm  and  deviation  <0.5 
mrad  at  attenuations  >35  dB. 

We  designed  the  attenuator  to  meet  a  number  of  cw 
and  pulsed  laser  attenuation  requirements  in  our  lab¬ 
oratory.  We  expect  that  the  wedged-plate  attenuator 
has  many  uses,  including  nonlinear  optical  studies, 
detector  linearity  and  calibration  measurements,  ab¬ 
solute  phototube  calibration,  tunable  laser  intensity 
control,  and  in  absorption  spectroscopy  as  part  of  a  null 
detector  servo  system. 
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We  have  invented  a  quasi-Gaussian  profile-transmittance  filter  based  on  radially  varying  the  phase  retardation 
in  a  birefringent  element.  The  radial  birefringent  element  has  been  applied  to  resonator  design  and  has  demon¬ 
strated  its  usefulness  in  generating  an  improved  resonator  spatial-mode  profile. 


It  has  been  known  for  more  than  one  decade  that  reso¬ 
nators  constructed  with  a  Gaussian  profile-reflector 
element  offer  the  possibility  of  arbitrary  large  mode 
diameters,*'4  thus  overcoming  the  small-mode-volume 
limitation  of  TEMoo  mode  resonators.  To  date,  ex¬ 
perimental  verification  of  the  theory  has  not  been 
demonstrated.  It  has  also  been  recognized  that  un¬ 
stable  resonators5  offer  an  advantage  of  large  mode 
volume,  but  at  a  sacrifice  of  mode  quality  because  of 
aperture-generated  Fresnel  fringes.6-7  Furthermore, 
the  stability  of  unstable  resonators  improves  with  in¬ 
creasing  magnification  and  thus  with  higher  output 
coupling.  Thus  unstable  resonators  have  been  limited 
to  relatively  high-gain  laser  media,  such  as  C02  (Ref.  8) 
and  Nd:YAG.9 

We  propose  an  elegant  yet  simple  optical  device,  the 
radial  birefringent  element  (RBE),  as  a  method  of 
realizing  a  Gaussian  reflectance  profile.  We  show  that 
the  RBE  can  be  used  in  both  stable  and  unstable  reso¬ 
nators  to  realize  the  advantages  predicted  by  theory.  In 
particular,  the  RBE,  in  conjunction  with  an  unstable 
resonator,  provides  a  smooth  spatial-mode  profile  and 
arbitrary  output  coupling  and  yet  maintains  the  sta¬ 
bility  associated  with  large-magnification  unstable 
resonators. 

The  RBE  is  shown  schematically  in  Fig.  1  as  the  end 
reflector  of  a  resonator.  The  element  consists  of  a  bi¬ 
refringent  crystal  with  a  radius  of  curvature  polished 
on  one  surface  situated  between  a  polarizer  and  an  end 
mirror.  The  reflectance  at  the  center  of  the  element  is 
determined  by  the  net  phase  retardation  <f>o  - 
2srAn/0A,  where  An  is  the  birefringence,  lo  is  the 
physical  thickness,  and  A  is  the  wavelength.  Since  the 
phase  retardation  #(r)  is  a  function  of  radial  position, 
the  reflectance  also  varies  radially.  If  we  assume  that 
the  reflected  beam  retraces  the  optical  path  of  the  for¬ 
ward  beam  at  the  element,  an  analysis  of  the  mirror, 
polarizer,  and  RBE  combination  using  Jones’s  matrix 
approach"*  shows  that  the  reflectance  and  transmit¬ 
tance  of  the  device  are  given  by 

R  ■  cos’-  0(r)  +  sin*  0(r)cos2  20,  (la) 

T  ■  sin2  ^(r)sin2  20,  (lb) 

when*  0  is  the  angle  between  the  polarizer  and  the 


principal  axis  of  the  birefringent  plate.  To  a  good  ap¬ 
proximation, 


2wAn 

A 


Hr) 


(2) 


where  r  is  the  radial  position  and  p  is  the  curvature,  with 
the  upper  and  lower  signs  corresponding  to  convex  or 
concave  curvatures,  respectively.  Equation  (la)  can  be 
written  in  the  form 


R(r )  =  cos2  20  +  sin2  20  cos2  <t>(r). 

which  explicitly  shows  that  the  reflectance  is  given  by 
a  constant  factor  depending  on  0  modulated  by  a  term 
varying  as  cos2  <t>(r)  with  amplitude  sin2  20.  For  0  =  45°, 
we  find  that  R{r )  =  cos2  <t>(r),  or 

prM-  ,3) 

If  we  want  to  set  R(r)  =  0  at  a  specific  radius  r(>,  de¬ 
termined,  for  example,  by  an  aperture  or  rod  radius 
within  the  resonator,  then  Eq.  (3)  can  be  inverted  to 
give 

lo  =  2xAn  +  mxl> 

p  =  r°2^2f  |cos-*|R(r0)]'/2  -  cos-*[R(0)],/2|-\ 

(4) 

where  m  is  an  integer. 
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Fig.  1.  The  radial  birefringent  element  consisting  of  a  po¬ 
larizer  and  a  birefringent  plate  with  curvature  p,  within  an 
unstable  resonator  cavity. 
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Fig.  2.  Reflectance  profiles  of  the  RBE  using  two  crystal  - 
quartz  elements.  Two  elements  have  equal  but  opposite 
curvature  p.  0j  and  are  the  angles  between  the  principal 
axis  of  the  quartz  elements  and  the  polarizer  axis.  The  mode 
radius  r0  is  chosen  equal  to  the  radius  of  the  Nd.YAG  rod. 


It  is  convenient  to  define  the  equivalent  Gaussian 
spot  size  of  the  RBE  wm  as 


,  _  ,  cos~1[/?(0)/e]1/2  -  cos~1lfi(0)J1/2 

0  cos-«[R(r0)J1/2-cos-Mfl(0)J^’  (  ) 


where  e  is  the  Neper  number.  Note  that  for  /?(0)  =  1 
and  R(ro)  =  0,  wm2  =  0.58  w02,  and  that,  unlike  for  a 
Gaussian  reflector,  weean  arrange  R(r0)  >  R( 0).  Thus 
the  RBE  is  a  pseudo-Gaussian  mirror  with  reflec¬ 


tance 


y/Rl o)c-<r*W) 

Equation  (5)  can  be  used  with  previous  analyses  of 
resonators1  with  Gaussian  reflectance  elements  as  the 
effective  spot  size  of  the  Gaussian  reflector. 

Finally,  for  ABCD  resonator  analysis  of  the  cavity 
shown  in  Fig.  1,  the  radius  of  curvature  of  the  equivalent 
mirror  is  given  by 


where  f  is  the  focal  length  of  the  birefringent  plate  with 
curvature  p  and  Rm  is  the  end  mirror’s  radius  of  cur¬ 
vature.  This  is  a  good  approximation  when  the  element 
is  a  distance  from  the  end  mirror  that  is  small  compared 
with  the  mirror’s  curvature  or  the  element’s  focal 
length. 

We  have  also  designed,  constructed,  and  analyzed11 
other  variations  of  the  RBE.  Of  interest  for  resonator 
studies  is  a  two-element  RBE  consisting  of  a  positive- 
curvature  element  next  to  a  negative-curvature  element 
By  varying  the  angles  0]  and  02,  which  are  the  angles 
between  the  polarizer  and  the  principal  axis  of  the  bi¬ 
refringent  elements,  various  reflectance  profiles  can  be 
generated.  Figures  2  shows  some  reflectance  profiles 
for  various  02  and  fixed  0j.  Here  both  elements  have 
equal  refractive  indices  and  equal  but  opposite  curva¬ 


tures  p.  This  variation  of  the  RRE  was  used  in  the 
resonator  studies  described  next. 

We  have  constructed  both  stable  and  unstable  reso¬ 
nators  using  the  RBE  for  a  Nd.  YAG  gain  medium.  The 
resonators  with  the  geometric  magnification  larger  than 
2  demonstrated  that  large  TEMoo  mode  spot  sizes  could 
be  supported,  as  expected  from  previous  theoretical 
analysis.1-3’4  However,  when  the  mode  size  is  big 
enough  to  extract  high  energy  from  the  rod,  the  dif¬ 
fraction  loss  that  is  due  to  the  rod  causes  the  stability 
of  such  a  resonator  to  be  marginal. 

The  RBE  in  combination  with  an  unstable  resonator 
offers  the  advantages  of  good  alignment  stability,  ex¬ 
cellent  output  spatial-mode  profile,  and  good  discrim¬ 
ination  against  higher-order  modes.  The  RBE  was 
constructed  by  combining  a  convex  and  a  concave 
crystal  quartz  lens  with  equal  p  =  25-cm  curvatures. 
The  mirror  curvatures  of  the  RBE  resonator  were  cho¬ 
sen  to  form  an  unstable  resonator  with  the  thermal  fo 
cusing  of  the  rod  taken  into  account.  Best  collimation 
at  the  output  was  achieved  with  the  combination  of  a 
-400-cm-curvature  mirror  and  a  — 200-cm-curvature 
mirror.  The  resonator  length  was  65  cm.  In  this  res¬ 
onator,  the  output  coupling  was  changed  by  rotating  the 
RBE.  The  optimum  performance  was  found  when  6\ 

=  45°  and  02  =  100°,  although  the  angles  were  relatively 
uncritical  with  A0j  =  ±5°  and  A02  =  dkl0°. 

We  chose  to  compare  the  RBE  resonator  performance 
against  the  standard  60-cm-long,  M  -  3.3  Nd:YAG 
confocal  unstable  resonator.9  Figure  3  shows  the 
transverse-mode  quality  of  the  RBE  unstable  resonator 
and  the  confocal  unstable  resonator  with  a  1.8-mm- 
diameter  polka-dot-coated  output  coupler.  The  video 
scans  clearly  show  the  flat-top  profile  of  the  RBE  res¬ 
onator  compared  with  strongly  Fresnel-modulated 
unstable-resonator  output.7  The  RBE,  with  its 
cos2  reflectance  profile,  effectively  eliminates 


RBE  UNSTABLE  RESONATOR 


CONFOCAL  UNSTABLE  RESONATOR 


Fig.  3.  Transverse-mode  profiles  of  the  RBE  and  standard 
confocal  unstable  resonator  configurations.  The  Polaroid 
exposures  of  the  output  are  at  the  right  side  and  the  vidicon- 
scanned  profiles  are  at  the  left  side.  Hie  slightly  asymmetric 
patterns  are  due  to  a  damage  spot  on  one  end  of  the  Nd:YAG 
rod. 
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Fig.  4.  Second-harmonic-generation  conversion  efficiency 
versus  input  pump  energy  for  the  confocal  unstable  resonator 
and  for  the  RBE  unstable  resonator  in  Type  II  KD*P. 


Fresnel  fringes  that  arise  from  diffraction  caused  by 
hard  apertures.  The  weak  Fresnel  fringes  on  the  profile 
of  the  RBE  resonator  are  due  to  the  finite  aperture  size 
of  the  YAG  rod. 

As  a  further  comparison  of  performance,  we  carried 
out  second-harmonic-generation  (SHG)  conversion- 
efficiency  measurements  in  a  Type  II  angle-phase- 
matched  KD*P  crystal.  Figure  4  shows  the  improved 
SHG  performance  of  the  RBE  unstable  resonator. 

In  conclusion,  we  have  invented  a  new  class  of  optical 
elements  based  on  radially  varying  the  phase  retarda¬ 
tion  in  a  birefringent  plate.  In  particular,  the  radial 
birefringent  element  has  been  applied  to  resonator  de¬ 
sign  and  has  demonstrated  its  usefulness  in  generating 
an  improved  resonator  spatial-mode  profile  and  im¬ 
proved  nonlinear  conversion.  The  RBE  is  simple  to 
implement  and  should  prove  useful  in  designing  un¬ 
stable  resonators  for  large-aperture  but  low-gain  laser 
media.  The  RBE  should  also  find  applications  as  a 
spatial  filter  or  soft  aperture. 
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ABSTRACT 

A  reliable  single  ended  remote  monitoring  LIDAR  operating  in  the  infrared  is 
described.  The  system  uses  a  1.4  to  4.0  pm  tunable  0P0  source  as  transmitter.  Con¬ 
tinuous  remote  monitoring  of  atmospheric  methane  over  a  17  hour  period  demonstrates 
the  reliability  of  the  system.  Simultaneous  remote  measurements  of  temperature  and 
humidity  using  absorption  lines  of  the  1.9  pm  water  band  achieve  relative  accuracies 
of  1 . 0°C  and  1  V.  respectively  over  a  45  sec.  averaging  period.  The  expected  sensitiv¬ 
ity  for  the  measurement  of  other  pollutants  with  absorption  lines  uithin  the  tuning 


range  of  the  transmitter  is  discussed. 


1.  INTRODUCTION 

The  need  for  single  ended  air  pollution 
monitoring  systems  is  now  well  estab¬ 
lished  and  the  methodology  of  making 
remote  measurements  is  understood1.  The 
monitoring  of  industrial  pollution  and 
the  measurement  of  trace  constituents 
are  efficiently  done  using  a  single- 
ended  remote  probing  technique.  Fur- 
thermore*  remote  probing  is  evolving 
into  an  important  research  tool  for 
studying  transport,  mixing  and  chemistry 
of  trace  constituents  in  the  atmosphere. 
Recent  experiments  have  also  shoun  the 
potential  for  simultaneous  remote  humid¬ 
ity  and  temperature  measurements  for 
metereology. 

In  their  review  article  K i 1  dal  and  Byer2 
found  that  remote  measurements  using 
absorption  promised  the  highest  sensi¬ 
tivity  and  lowest  required  powerlevels. 
In  a  later  paper.  Oyer  and  Garbuny3  ana¬ 
lyzed  remote  measurement  methods  includ¬ 
ing  single  ended  absorption  with  topo¬ 
graphic  targets  and  depthresol ved 
absorption  using  Mie  scattering  as  a 
distributed  retroref lector.  For  long 


path  absorption  a  narrow  bandwidth  laser 
beam  is  transmitted  through  the  atmos¬ 
phere  to  a  distant  topographic  target. 
Some  portion  p  of  this  beam  is  scattered 
back  by  the  non  cooperative  target,  and 
is  collected  by  a  telescope  at  the 
transmitter  end.  The  received  energy 
SCv)  is  determined  by  using  conservation 
of  energy  leading  to  the  LIDAR-equation 

R 

-2<f{ascaCr)+NCr)o(v)}dr 
SCv)=E0Cv)flA/R2p*e  0 

Ml 

where 

Eo=transmi tted  light-energy  of  fre¬ 
quency  v 

inefficiency  of  optical  system 
A=area  of  receiving  telescope 
R=distance  to  backscattering  target 
p(R)=backscattering  due  to  target  at 
distance  R 

«kca=attenuation  coefficient  due  to 
atmospheric  scattering-losses 
N(r)=numberdensity  of  species  to  be 
probed 

g(v)=absorption  cross  section  of 
species  at  wavelength  v 
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In  the  near  infrared*  topographic  tar¬ 
gets  scatter  1  to  10%  of  the  incident 
light  back  toward  the  receiving  tele¬ 
scope.  The  average  concentration  of  a 
pollutant  along  the  lightpath 
Ns/M(  r)dr/R  is  measured  by  transmitting 
two  pulses  with  different  frequencies  v. 
If  the  frequency  change  is  small*  only 
the  absorption  cross  section  of  the  pol¬ 
lutant  changes*  so  that  the  concentra¬ 
tion  N  derived  from  (1)  is 


ln{S(vi)/S(v2)J 

N=  - 

2R{«(V2)-c(Vi)} 


12) 


where  the  absorption  cross  sections 
o(Vi)  have  been  previously  determined. 
Fixed  frequency  lasers  have  only  a  few 
coincidences  with  absorption  lines  of 
interest'  so  that  for  most  applications 
a  tunable  laser  is  preferred. 

From  1974  on  many  remote  monitoring  sys¬ 
tems  have  been  developed  and  tested 
using  tunable  dye  lasers5'6.  Initially 
the  reliability  of  many  of  these  systems 
was  marginal,  but  recently  a  number  of 
systems  with  excellent  reliability  using 
commercially  available  Nd:Yag  pumped  dye 
lasers  have  been  reported7'8.  A  series 
of  field  measurements  have  demonstrated 
the  usefulness  of  remote  pollution 
detection8.  However,  dye  laser  measure¬ 
ments  in  the  visible  and  ultra-violet 
region  is  limited  to  the  measurement  of 
only  three  pollutants:  SO2,  NO2,  and  O3. 


To  extend  the  LIDAR  measurements  to 
other  atmospheric  constituents,  the 
infrared  spectral  region  must  be  used. 
Line  tunable  C02*  recently  doubled  CO2 
and  OF  lasers  have  been  used8'’0'’’,  but 
are  limited  to  spectral  coincidences 
uith  absorption  lines  of  pollutant  mole¬ 
cules.  To  achieve  wider  pollutant 
selection  and  better  measurement  sensi¬ 
tivity*  we  have  developed  a  continuously 
tunable  high  power  infrared  source  which 
is  used  as  transmitter  for  a  single 
ended  remote  measurement  system.  The 
transmitter  is  an  Optical  Parametric 
Oscillator  (0T0)»  which  is  continuously 


tunable  from  1.4  to  4.0  pa.  Initial 
measurements  of  S02*2  and  also  CH*’3 
with  this  system  have  been  reported  by 
Baumgartner  and  Byer.  Improvements  in 
the  design  of  the  0P0  transmitter  which 
increased  the  reliability  and  spectral 
resolution  considerably  were  described 
by  Brosnan  and  Dyer1'.  Increased  output 
energy  up  to  70  mJ  per  pulse  at  10  Hz 
was  achieved  uith  an  Optical  Parametric 
Amplifier  (OPA)  discussed  by  Baumgartner 
and  Byer15.  Extensive  software  devel- 
opement  has  enabled  automatic  pollution 
measurements  to  be  made  over  timescales 
of  many  hours  as  demonstrated  by  recent 
CH*  measurements’6. 


II.  SYSTEM  DESCRIPTION 

The  Stanford  LIDAR  transmitter  is  an 
angle  phasematched,  singly  resonant 
LiNb03  Optical  Parametric  Oscillator 
(OPO).  For  the  methane  measurements  the 
0P0  was  amplified  by  an  Optical  Parame¬ 
tric  Amplifier  (OPA).  OPO  and  OPA  are 
described  in  detail  by  Brosnan  and 
Byer"  and  Baumgartner  and  Byer’5. 

Fig.l  shows  a  diagram  of  the  transmit¬ 
ter.  The  Q-suitched  unstable  resonator 
Nd:Yag  laser  generates  220  mJ  energy  per 
pulse  with  a  10  Hz  repetition  rate.  The 
length  of  the  resonator  is  1.40  m  to 
generate  18  nsec  long  Q-suitched  pulses. 
This  pulselength  is  ideal  to  pump  the 
OPO,  since  the  optical  intensity  of  the 
pumpbeam  is  low  enough  to  avoid  damage 
of  the  LiHbOs  crystal  due  to  self  focus¬ 
ing,  but  the  oscil lotion  threshold  of 
the  OPO  still  is  lou  due  to  decreasing 
build-up  losses". 

The  output  of  the  Hd:Yag  laser  is 
slightly  riemagnified  to  a  5  mm  diameter 
beam  and  directed  into  a  folded  15  m 
long  path.  This  distance  is  adequate* 
to  transform  the  unstable  resonator  1.06 
pm  beam  into  the  far  field  Airy-disk 
made  with  CO#  of  the  beam  energy  in  the 
centra]  lobe.  Oscillation  threshold  of 
the  OPO  is  considerably  decreased  with 
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this  pump  beam  compared  to  pumping  uith 
the  near  field  unstable  mode.  The  dan¬ 
ger  of  surface  damage  to  the  LiNbOs 
crystal  is  also  removed.  The  pointing 
stability  of  the  long-cavity  MdsYag 
laser  and  the  far  field  converter  is 
critical  for  the  reliability  of  the 
system.  Me  used  fused  silica  rods  to 
stabilize  the  Nd:Yag  resonator  struc¬ 
ture#  and  mounted  the  steering  mirrors 
for  the  folded  15  m  path  length  far 
field  converter  directly  onto  the  laser 
resonator  structure.  This  arrangement 
gives  excellent  pointing  stablilty  of 
the  beam. 

For  operation  of  the  OPA  40 '/,  of  the  130 
mo  of  1.06  pm  energy  is  split  off  into  a 
Nd:Yag  amplifier,  uhile  the  mainpart  is 
directed  through  a  2.2X  demagnifying 
telescope  to  pump  the  0P0.  The  0P0  uses 
a  50  mm  long#  20  mm  diameter  LiNb03 
crystal  to  generate  the  signal  and  idler 
output.  The  signal  wave  tunes  from  2.1 
to  1.4  pm,  uhile  the  idler  wave  tunes 
from  2.1  to  4.0  pm  for  a  4°  change  in 
the  LillbOa  phase  matching  angle.  The 
tuning  range  is  limited  by  absorption  of 
the  idler  wave  in  the  LiNb03  crystal  at 
4  pm.  The  0P0  uses  an  L-shaped  resona¬ 
tor  uhich  reflects  only  the  signal  wave. 
The  resonator  uses  a  flat  outputcoupler 
uith  502  reflectivity  between  1.4  and 
2.1  pm,  .  and  a  600  gr/mm  grating.  In 
order  to  improve  the  resolution  of  the 
grating,  a  10:1  ZnSe  prism  beam 
expander17  is  inserted  into  the  cavity. 
With  grating  and  prism  beam  expander  the 
lineuidth  at  half  intensity  points  is 
less  than  0.6  cm'1  for  the  signal  uave. 
To  improve  the  resolution  further,  a 
tilted  etalon  can  be  inserted  into  the 
cavity.  For  a  2  mm  thick  etalon  the 
signal  output  lineuidth  is  less  than  0.1 
cm'1*  The  lineuidth  of  the  idler  output 
Is  about  equal  to  the  sum  of  the  widths 
of  signal  and  pump  beam.  The  Nd:Yag 
laser  has  a  lineuidth  of  0.4  cm*1,  but 
can  easily  be  reduced  to  0.01  cm'1  by 
inserting  a  tilted  etalon. 

Output  energy  of  the  0P0  is  5  mJ  in  the 
signal  wave  over  most  of  the  tuning 


range.  Due  to  the  absorption  losses  of 
the  idler  uave  in  the  LiKbOs  crystal, 
the  threshold  increases  near  the  end  of 
the  tuning  range  and  the  output  energy 
decreases.  The  energy  of  the  idler  out¬ 
put  is  equal  to  the  energy  in  the  signal 
uave  multiplied  by  the  ratio  of  the  pho¬ 
ton  energies. 

Astigmatic  beam  divergence  uas  observed 
in  the  0P0  output  beams  when  the  grat¬ 
ing-prism  beam  expander  is  used  for 
lineuidth  control.  With  the  etalon  the 
output  beam  has  a  symmetrical  divergence 
and  can  be  collimated. 

The  OPA  consists  of  a  Nd:Yag  pumped  50mm 
x  20  mm  angle  phasematched  LiNbOs  crys¬ 
tal.  It  increases  the  output  energy 
from  the  0P0  to  more  than  12  mJ  in  the 
signal  uave.  Houever,  for  the  tempera¬ 
ture  and  humidity  measurements  only  the 
0P0  uas  used.  This  simplified  the  set 
up,  and  yet  provided  enough  energy  to 
see  a  strong  backscatter ing  signal  from 
topographic  targets.  Table  1  summarizes 
the  transmitter  parameters  for  the  dif¬ 
ferent  measurements. 

About  IX  of  the  tunable  light  energy  is 
used  for  diagnostic  purposes.  The  lack 
of  wavelength  resettabi 1 i ty  of  the  OPO 
in  early  experiments  led  to  the  intro¬ 
duction  of  a  photo-acoustic  cell  to  mon¬ 
itor  the  proper  tuning  of  the  OPO  by 
observing  the  absorption  cross  section 
during  an  extended  measurement. 
Improvements  of  the  tuning  mechanism 
made  continuous  monitoring  of  the  cross 
section  superfluous  and  later  the  cell 
uas  used  only  for  spec-troscopic  measure¬ 
ments.  An  InSb  detector  is  used  as  ref¬ 
erence  detector  to  ratio  out  energy 
fluctuations  of  the  0P0/0PA.  An  absorp¬ 
tion  cell  follows  to  determine  the 
absorption  cross  section  of  a  pollutant 
at  a  given  wavelength. 

The  main  portion  of  the  beam  is  expanded 
in  a  10:1  telescope  to  about  2.5  cm 
diameter  and  transmitted  coaxially  from 
the  receiving  telescope  on  the  roof  of 
our  laboratory.  This  receiving  tele- 
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scope  is  a  40  cm  diameter,  f =3  Newtonian 
design.  It  collects  and  focuses  the 
backscattered  light  onto  a  liquid  nitro¬ 
gen  cooled  InSb  detector. 

Fig. 2  shous  a  diagram  of  the  electron¬ 
ics.  The  amplified  signal  and  reference 
pulses  are  integrated  and  digitized  for 
the  computer,  which  ratios  and  averages 
the  signals.  The  signal  from  the  photo¬ 
acoustic  cell  is  amplified  in  a  Lock-in 
amplifier,  which  is  tuned  to  the  reso¬ 
nance  frequency  of  the  cell  and  syn¬ 
chronized  from  the  trigger  signal  of  the 
laser. 

Uavelength  selection  of  the  0P0  is  con¬ 
trolled  by  a  PDP-11/10  minicomputer. 
Crystal  and  grating  are  set  according  to 
values  in  a  table,  while  the  position  of 
the  etalon  is  calculated  from  a  set  of 
parameters.  An:  interactive  program  was 
developed  to  simplify  operation  of  the 
OPO  LIDAR,  and  to  make  automatic  meas¬ 
urements  possible.  Besides  service  rou¬ 
tines  for  calibration  and  alignment  of 
the  OPO,  the  program  enables  automatic 
scans  of  the  OPO,  increasing  or  decreas¬ 
ing  the  uavelength  in  predetermined 
steps.  In  another  mode  the  program 
tunes  the  OPO  between  two  or  three  pre¬ 
determined  wavelengths,  averages  the 
received  signals  from  each  uavelenth  and 
also  calculates  the  rms-variation  of 
this  average  for  later  error  analysis. 
It  then  calculates  the  concentration  of 
a  pollutant  using  eqn.(2).  This  value 
is  displayed  in  real  time  as  well  as 
stored  for  later  use.  This  control  rou¬ 
tine  proved  to  be  essential  for  the 
measurements  reported  in  later  sections 
of  this  paper. 

Without  the  etalon  in  the  cavity,  the 
OPO  can  be  scanned  under  computer  con¬ 
trol  continuously  from  1.4  to  4.0  pm,  a 
range  of  more  than  4500  cm'*.  Fig. 3 
shous  an  example  of  the  scan  capabili¬ 
ties  of  the  OPO.  In  the  lower  trace  it 
shous  an  atmospheric  transmission  spec¬ 
trum  from  31  May  1979.  We  used  a  build¬ 
ing  775  m  away  as  backscntter ing  target, 
so  the  total  atmospheric  pathlength  is 


1.55  km.  The  scan  covers  a  450  cm'1 
range  from  2850  to  3300  cm*1,  which 
includes  numerous  absorption  lines 
mainly  due  to  water  vapor  and  methane, 
although  many  other  trace  constituents 
have  absorption  lines  in  this  region. 
To  confirm  the  accuracy  of  our  scan  and 
to  identify  lines  we  use  the  AFCRL  data 
base1*.  The  upper  trace  in  fig. 3  shous 
a  simulated  transmission  spectrum  using 
the  AFCRL  data.  For  this  simulation  ue 
assumed  6  torr  II2O  and  3  ppm  CH%  in  the 
atmosphere.  The  lineuidth  was  assumed 
to  be  1.5  cm"1  However,  the  program  that 
calculates  the  transmission  from  the 
AFCRL  data  uses  a  rectangular  slit  func¬ 
tion,  which  is  not  a  very  good  approxi¬ 
mation  of  the  OPO  lineshape,  and  results 
in  discrepancies  betueen  the  calculated 
and  the  observed  spectrum.  But  the  gen¬ 
eral  agreement  of  both  traces  in  fig. 3 
confirms  the  usefulness  of  the  OPO- 
source  for  spectroscopic  applications. 

It  is  considerable  more  complex  to  con¬ 
trol  continuous  scans  with  the  etalon 
inserted  into  the  OPO  cavity.  Grating 
and  etalon  position  have  to  be  synchron¬ 
ized  to  within  0.1  cm*1.  This  requires 
excellent  stability  of  the  tuning  mecha¬ 
nism.  Furthermore,  if  a  scan  over  more 
than  one  free  spectral  range  of  the 
etalon  is  required,  the  etalon  has  to  be 
reset  to  its  zero  position  during  the 
scan.  The  length  of  a  scan  is  deter¬ 
mined  by  the  accuracy  with  which  the 
etalon  free  spectral  range  and  the  tun¬ 
ing  rate  of  the  grating  are  known.  Me 
developed  extensive  software  to  measure 
semi-automatical ly  the  etalon  parameters 
and  so  enable  scans  .aver  many  etalon 
free  spectral  ranges.  We  demonstrated 
scans  over  more  than  10  free  spectral 
ranges,  which  in  the  case  of  a  2  mm 
thick  fused  silica  etalon  is  about  20 
cm" 1 .  For  our  applications,  this  scan 
range  is  more  than  adequate.  Tig. 4  dem¬ 
onstrates  the  scanning  capabilities  and 
the  linewidth  of  our  OPO  source  with  and 
without  intracavity  etalon.  It  shows 
the  spectrum  of  the  Q-branch  of  the  2va 
overtone  transition  in  CIl*  at  lou  pres¬ 
sures  obtained  in  a  photo  acoustic  cell 
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(see  also  the  high  resolution  spectrum 
obtained  by  Bobin"l.  The  resolution  of 
the  0P0  uith  and  without  the  etalon 
inserted  is  apparent.  The  etalon  nar¬ 
rowed  the  0P0  linewidth  to  0.08  cm*1. 

The  reliability  of  the  Nd:Yag  OPO  trans¬ 
mitter  has  been  demonstrated  by  more 
than  a  year  of  operation  without  damage 
to  the  pumplaser  or  the  OPO.  Since  the 
LIDAR  consists  only  of  solid  state  com¬ 
ponents.  no  periodic  maintenance  is 
required. 


[THANE  MEASUREMENT; 


Me  selected  the  measurement  of  methane 
concentration  to  demonstrate  the  relia¬ 
bility  of  our  system.  The  atmospheric 
methane  concentration  is  recorded  at  the 
Bay  Area  Air  Quality  Management  District 
point  monitoring  station  in  Redwood  City 
10  km  from  our  laboratory.  This  allows 
the  possibility  of  comparing  the  accu¬ 
racy  of  our  system  with  an  independent 
station. 


Fig. 5  shows  two  transmission  spectra 
near  the  CH*  absorption  region  at  3.4 
pm.  Both  spectra  were  obtained  with  the 
idler  output  of  the  OPO/OPA.  The  lineu- 
idth  of  the  output  was  1.4  cm*'.  The 
upper  trace  shows  the  atmospheric  trans¬ 
mission  betueen  our  system  and  a  build¬ 
ing  at  775  m  distance,  which  was  used  as 
non  cooperative  backscattering  target. 
Both  CH*  and  H2O  absorption  lines  are 
present.  To  clearly  identify  the  CHi, 
lines  and  to  aid  in  selecting  a  trans¬ 
ition  that  is  free  of  water  vapor  inter¬ 
ference.  an  absorption  spectrum  of  CH* 
was  measured  in  the  laboratory.  A  sec¬ 
tion.  of  the  scan  centered  on  the 
P-branch  is  also  shown  in  the  lower 
trace  in  Fig. 5.  The  two  scans  in  this 
figure  can  be  overlaid  to  show  that  the 
P<2),  P<7) ,  P  C  9 )  and  PC  10>  methane 
transitions  are  nearly  water  vapor  free. 
Additional  confirmation  is  also  obtained 
from  the  spectroscopic  information  in 
the  AFCRL  data".  The  ability  to  select 


interference-free  absorption  bands  and 
to  verify  -the  location  of  spectral  lines 
by  scanning  the  spectrum  is  a  major 
advantage  of  the  tunable  source. 

Me  selected  the  P(10>  transition  for  our 
measurements  because  it  has  the  proper 
absorption  strength  for  CH*  measurements 
over  the  selected  775  m  long  path.  For 
measurements  over  shorter  paths.  the 
P(7)  line  with  a  much  larger  cross  sec¬ 
tion  can  be  used.  The  computer  was  pro¬ 
grammed  to  tune  on  and  off  the  PC  103 
line  every  20  pulses  or  2  seconds.  The 
data  were  normalized  and  processed  to 
show  the  concentration  of  methane.  An 
effective  cross  section  for  the  CH« 
P C 1 0 3  line  was  used  based  on  the  resolu¬ 
tion  of  the  OPO/OPA  source.  To  correct 
for  any  long-term  drift,  the  cross  sec¬ 
tion  was  measured  continuously  during  a 
run  by  the  photo  acoustic  cell.  The 
measured  value  of  the  cross  section  was 
used  to  calculate  the  methane  concentra¬ 
tion  in  real  time. 

Fig. 6  shows  a  segment  of  a  17  h  CH« 
measurement  run.  Here,  the  raw  measure¬ 
ments  are  shoun  for  20  shots/line  aver¬ 
aging  or  4  sec  per  point.  The  rms  fluc¬ 
tuations  are  calculated  to  be  0.6  ppm. 
To  improve  the  precision  of  the  measure¬ 
ment,  longer  averaging  times  can  be 
used.  The  heavy  line  in  Fig. 6  shows  the 
same  data  after  later  averaging  to 
increase  the  relative  accuracy  to  0.09 
ppm.  The  time  resolution  is  80  sec.  per 
point  in  this  case. 

Fig. 7  shows  the  full  17  hour  long  meas¬ 
urement  during  which  the  system  was 
operating  automatically  under  computer 
control.  During  the  measurements  con¬ 
cern  developed  that  the  mesasured  CHi, 
levels  were  high  since  the  global  avei — 
age  CH*  level  is  close  to  1.5  ppm. 
Hourly  CHi,  average  measurements  were 
obtained  from  the  Redwood  City  Bay  Area 
Air  Quality  Management  District  point 
monitoring  station  10  km  from  our  lidar 
at  Stanford.  The  concentration  measure¬ 
ments  from  this  point  monitoring  station 
are  indicated  by  open  circles  in  Fig. 7. 
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The  agreement  between  the  concentration- 
values  from  both  measurements  is  as  good 
as  ft  .can  be  expected  considering  the 
physical  seperation  between  Stanford 
Lidar  system  and  the  Redwood  City  moni¬ 
toring  station. 

To  check  further  on  the  correlation 
betueen  the  Stanford  and  the  Redwood 
City  station  measurements.  a  second 
extended  CH*  measurement  was  made  on 
June  6.  1979.  Fig. 8  shows  the  results. 
In  this  case,  the  wind  was  blowing  down 
the  bay  from  Redwood  City  toward  Stan¬ 
ford.  The  wind  velocity  and  distance 
gave  a  20-  to  30-  min  delay  time  from 
Redwood  City  to  Stanford.  This  delay 
correlates  with  the  peak  CH*  mesurement 
at  2200  at  Redwood  City  and  at  Stanford 
at  2230.  Later  in  the  evening  when  the 
wind  ceased.  Stanford  measurements 
became  independent  of  the  Redwood  City 
station  measurements. 

These  measurements  demonstrate  the  reli¬ 
ability  and  accuracy  of  remote  air  pol¬ 
lution  measurements.  Other  pollutants 
with  absorption  lines  within  the  tuning 
range  of  the  0P0  source  can  be  measured 
with  the  same  relative  accuracy.  The 
measurement  of  other  species  is  dis¬ 
cussed  in  section  V. 


IV.  TEMPERATURE  MEASUREMENT 

Tor  accurate  pollutant  measurements  it 
is  important  to  determine  the  tempera¬ 
ture  of  the  pollutant,  since  the  absorp¬ 
tion  cross  section  is  temperature 
dependent.  This  effect  is  usually  neg¬ 
lected  for  the  measurement  of  pollutants 
at  atmospheric  temperatures.  However 
for  emissions  from  sources  at  elevated 
temperatures  the  temperature  correction 
can  become  very  important.  By  properly 
selecting  appropriate  lines  of  a  molecu¬ 
lar  species  it  is  possible  to  make 
remote  temperature  and  concentration 
measurements  of  the  species  simultane¬ 
ously. 

A  number  of  groups  have  proposed  and 


analyzed  methods  of  making  remote  temp¬ 
erature  measurements,  and  preliminary 
remote  temperature  measurements  have 
been  reported  using  Raman  LIDAR20  and  a 
two  frequency  absorption  measurement21. 
Remote  temperature  measurements  using 
atmospheric  transmittance  measurements 
of  three  wavelengths  have  been  analyzed 
by  Mason22,  Schwemmer  and  Wilkerson23. 

The  temperature  dependence  of  the 
absorption  cross  section  is  given  by  the 
Boltzmann  distribution  and  a  partition 
function  Q(T)  as 

c(v?,T)=  c(v5,To)Q(T)exp{(E0i/k)Cl/T0-1/T)} 

[3] 

for  an  arbitrary  reference  temperature 
To,  where  E°i  is  the  lower  state  energy 
of  the  transition  at  frequency  v;,  and  k 
is  the  Boltzmann  constant. 

Measuring  the  absorbance 

A(Vi,T)  =N  a(vi,T)R 

for  two  lines  of  one  species  with  dif¬ 
ferent  lower  state  energy  E°i  ue  can  use 
eqn . ( 3)  to  determine  the  temperature  T  ! 


To 

T=  - 

kT0  A(vn,T)  ff(vi.T) 

1 - {ln( - )-ln( - )} 

E°i-E°2  Atvj , T)  o(v2,T0) 

Ml 

This  expression  simplifies  for  T-T&  <<To 
to 

kT02  ACvi.T) 

T= - *ln( - )+D 

E°,-E02  A(v2,T)  : 

(51 

with 

kTo2  cr(v2.T0) 

D= - «{,n( - )}+T0 

E°i-E°2  c(v1tTo) 

The  constant  D  can  be  found  cither  using 
spectroscopic  data,  or  it  can  be  deter¬ 
mined  empirically  with  calibrated  temp¬ 
erature  measurements.  Once  the 


temperature  and  so  the  cross  section  is 
known*  it  is  straightforward  to  calcu¬ 
late  the  density  of  the  species  using 
eqn.(2). 

The  accuracy  of  a  remote  temperature 
measurement  is  determined  by  the  signal 
to  noise  ratio  (SNR)  and  the  number  of 
averaged  signals  (n)  .  From  eqn.(5)  one 
finds  that 

kT0*  1  2  1+exp(2A(vi.T))  I 

ST  = - * - {  X  ( - )) 

E#1-E°2  SNR)/n  i  =  1  A2(vi*T) 

(6) 

ST  is  minimized  for  A(v;,T)  «1.1,  where 
Cqn.(6)  simplifies  to 


E°1-E°2  SNRi/n 
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For  an  atmospheric  temperature  measure¬ 
ment*  we  selected  H2O  as  absorbing  spec¬ 
ies  for  the  following  reasons: 

H2O  allows  simultaneous  measurement 
of  temperature  and  humidity,  an 
important  capacity  for  remote  meter- 
eological  measurements. 

Due  to  the  complexity  of  the  H2O 
spectrum,  nearby  lines  with  a  large 
lower  state  energy  difference  can  be 
found.  This  simplifies  the  measure¬ 
ment  since  the  transmitter  has  to  be 
tuned  only  over  small  wavelength 
increments. 

Lines  generating  the  appropriate 
absorbance  to  minimize  eqn.(6)  for 
different  pathlengths  can  be  found. 

Me  used  the  AFCRL-tapes1*  to  select 
appropriate  H20  lines  for  temperature 
and  humidity  measurements  over  a  755  m 
long  path  (for  a  discussion  of  this  line 
selection  see  ref. 241.  Me  selected  two 
lines  of  the  1.9  pm  H2O  band  at  5650.41 
cm*’  and  5651.33  cm*1.  The  difference 
in  lower  state  energy  is  1456  cm*1. 
From  eqn.(7),  we  estimate  a  temperature 


accuracy  of  1°C  for  a  product  of 
SNR\/n=175.  For  a  SNR  of  25,  this  means 
a  50  shot  average*  which  promises  a  good 
time  resolution  for  this  measurement. 

The  lower  trace  in  Fig. 9  shows  an  atmos¬ 
pheric  transmission  scan  with  a  0.1  cm*1 
resolution  tuned  over  a  17  cm*1  range 
around  the  selected  lines.  The  upper 
trace  is  a  simulated  transmission  spec¬ 
trum*  using  AFCRL-tapes1* .  Me  used 
scans  similar  to  this  to  calibrate  the 
tuning  of  the  0P0  and  to  identify  the 
temperature  sensitive  lines. 

Once  the  wavelength  of  these  lines  is 
determined  and  the  0P0  aligned*  measure¬ 
ments  of  temperature  and  humidity  are 
done  under  computer  control.  The  return 
signal  at  each  frequency  is  averaged 
over  several  shots  before  tuni.ng  to 
another  wavelength.  To  keep  errors  due 
to  water  vapor  fluctuations  small,  it  is 
desirable  to  tune  the  0P0  after  each 
shot  and  average  the  temperature  data 
afterwards.  However,  tuning  between 
different  wavelengths  takes  up  to  0.5 
sec.  To  improve  the  time  resolution  of 
our  measurements  we  average  30  shots  per 
line  before  tuning  to  the  next  line. 
Mhen  all  three  required  wavelengths  are 
probed,  temperature  and  humidity  are 
calculated.  This  cycle  continues,  until 
a  total  number  of  150  shots  per  line  is 
reached.  Then  the  computer  calculates 
the  average  of  five  temperature  and 
humidity  data  points  as  well  as  the 
error  due  to  return  signal  fluctuations, 
stores  and  displays  those  values. 

*  a 

Fig. 10  shows  the  record  of  temperature 
and  humidity  from  the  morning  of  the 
11th  March  1980.  Each  point  represents 
an  average  of  150  shots/line,  or  a  time 
resolution  of  45  sec.  The  outside 
weather  was  cloudy  at  first*  with  rain 
starting  at  0740.  The  rain  increased 
until  it  was  so  strong  at  0810  that  the 
measurement  had  to  be  interrupted.  At 
0820  the  rain  stopped*  the  cloud  cover 
slowly  dissipated,  and  the  wind 
increased.  The  remotely  measured  temp¬ 
erature  is  compared  with  a  record  of  a 


thermograph  at  the  telescope.  The 
general  agreement  is  good.  Some  discre¬ 
pancies  of  the  tuo  measurements  are 
expected  since  the  remote  measurement 
averages  over  a  path  from  10  to  60  m 
above  the  ground  and  775  m  length,  uhile 
the  thermograph  is  located  on  the  roof 
of  a  building.  The  rms-deviation  for 
each  temperature  measurement  is  calcu¬ 
lated  from  the  return  signal  fluctua¬ 
tions.  It  varies  with  the  metereol ogi- 
cal  conditions:  in  the  early  morning, 
the  SNR  uas  measured  to  be  25.  The  var¬ 
iance  of  the  temperature  measurements 
then  is  1.1°C,  which  is  only  slightly 
worse  than  the  estimate  in  eqn.(7).  The 
SNR  decreases  later  due  to  the  rain,  and 
accordingly  the  rms-deviations  increase 
to  1.6°C.  After  the  rain  stopped,  the 
SNR  improved  again  to  about  25.  How¬ 
ever.  the  temperature  uncertainety  stays 
at  1.5°C  The  reason  for  this  increased 
error  can  be  found  in  the  stronger 
humidity  fluctuations  due  to  the  uind. 
Any  variation  in  absolute  humidity  dur¬ 
ing  the  10  sec.  measurement  cycle  for 
each  temperature  point  increases  the 
uncertainety  of  this  measurement. 
Shorter  measurement  cycles  will  improve 
the  accuracy  of  the  measurement. 

The  calibration  of  our  temperature  data 
is  done  by  minimizing  the  difference 
between  the  remote  temperature  data  and 
the  thermograph  data.  Me  have  collected 
more  than  10  hours  of  temperature  data, 
and  have  determined  a  value  for  the  con¬ 
stant  D  as  35°C.  uhich  is  19°C  less  than 
predicted  from  the  available  spectro¬ 
scopic  data. 

Me  have  demonstrated  the  feasibility  to 
measure  temperature  and  humidity  simul¬ 
taneously.  The  1.9pm  H2O  band  makes 
temperature  measurements  with  1.0°C 
accuracy  possible.  For  the  present 
measurements  of  atmospheric  temperature 
the  selected  lines  are  optimal.  How¬ 
ever,  for  temperature  measurements  at 
elevated  temperatures  such  ns  encoun¬ 
tered  in  stack  exhausts  a  different 
selection  of  lines  has  to  be  made. 


JL. _ MEASUREMENT  OF  OTHER  ATMOSPHERIC 

TRACE  COHSIIIVENTS 

In  this  section  we  discuss  the  possibil¬ 
ity  of  measuring  other  trace  constitu¬ 
ents  in  the  atmosphere  with  the  OPO 
HOAR.  Most  pollutants  have  absorption 
bands  in  the  1.4  to  4.0  pm  tuning  range 
of  the  OPO.  Houever,  some  of  the  lines 
are  too  weak  to  measure  atmospheric  con¬ 
centrations,  while  others  lie  in  regions 
of  strong  water  vapor  absorption,  uhich 
makes  measurements  over  long  paths 
impossible. 

To  estimate  the  sensitivity  of  a  meas¬ 
urement  of  a  certain  species  it  is 
assumed  that  a  IX  change  in  transmission 
can  be  measured.  For  the  reported  meas¬ 
urements  this  is  true  for  averages 
between  150  laser  pulses  per  line  for 
the  H2O  measurements,  to  2000  pulses  per 
line  for  the  CH*  measurements,  uhich 
means  averaging  times  betueen  30  and  400 
sec.  The  variation  of  the  averaging 
time  is  due  to  the  variation  of  the  out¬ 
put  energy  and  pulse  to  pulse  stablilty 
of  the  OPO  over  its  tuning  range,  but 
also  due  to  water  vapor  continuum 
absorption  around  3.4  pm,#.  _ 

Me  discuss  measurement  capabilities  of 
the  major  pollutants  as  identified  by 
Mright  et  a),  in  ref. 25.  . .  Spectroscopic 
information  uas  obtained  from  different 
sources,  the  AFCRL  data  base18  being  the 
most  important  one  since  it  contains  the 
major  atmospheric  constituents  and  is 
ideally  suited  to  check  for  interfer¬ 
ences  of  different  species.  An  overview 
over  the  measurement  capabilities  is 
summarized  in  table  2. 

H*0:  The  influence  of  water  vapor  on 
atmospheric  processes  and  measurements 
is  so  important,  that  it  must  be  dis¬ 
cussed  first.  The  measurement  of  humid¬ 
ity  is  easily  possible  in  the  infrared, 
as  was  demonstrated  in  section  IV.  A 
large  number  of  lines  exists,  uhich 
makes  it  possible  to  select  one  with  the 
ideal  cross  section  for  n  specific  meas¬ 
urement.  The  two  water  vapor  absorption 
bands  of  importance  for  measurements 
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with  the  OPO  transmitter  are  the  1.9  pun 
and  2.7  ptm  bands.  These  bands  absorb  so 
strongly*  that  measurements  over  long 
paths  in  the  regions  from  1.8  to  2.0  urn 
and  2.4  to  3.2  |im  are  not  possible.  The 
surrounding  regions  contain  lines  of 
ideal  strengths  for  atmospheric  humidity 
measurements.  The  large  number  of  uater 
vapor  lines  leads  to  interference  prob¬ 
lems  in  the  measurements  of  many  other 
pollutants.  Thus  it  is  necessary  to 
measure  the  absolute  humidity  first  in 
order  to  correct  absorption  measurements 
on  pollutant  lines  influenced  by  nearby 
N2O  lines. 

CO2:  On  a  global  scale,  carbondioxide  is 
important  to  monitor  due  to  its  effect 
on  the  thermal  balance  of  the 
atmosphere26 .  On  a  local  scale,  CO2  can 
he  used  as  a  tracer  for  emissions  from 
combustion  processes,  and  to  measure  the 
temperature  of  these  emissions.  The 
atmospheric  background  of  C02  is  330  ppm 
in  an  unpolluted  atmosphere,  and  ranges 
up  to  350  ppm  in  polluted  regions.  For 
remote  measurements  with  the  OPO  trans¬ 
mitter,  the  vt+2v2+V3  -combination  band 
at  4983  cm*1  is  ideal.  Using  the  AFCRL 
data  base'®  and  our  own  atmospheric 
transmission  spectra  ue  determined  that 
the  P(S),  P( 10)  and  P(36)  lines  are 
interference  free  from  uater  vapor 
absorption.  These  lines  are  well  suited 
to  measure  atmospheric  background  as 
uell  as  source  emissions  of  CO2.  Using 
the  570  cm'1  lower  level  energy  differ¬ 
ence  of  the  P(8)  and  the  P(3G)  lines* 
temperature  of  sources  can  be  determined 
with  a  5°C  accuracy. 

CO:  Carbon  monoxide  is  knoun  to  be  a 
health  hazard  in  high  concentrations. 
The  Clean  Air  Act  sets  the  CO  exposure 
limit  to  an  average  value  of  less  than 
35  ppm  over  one  hour  or  9  ppm  over  24 
hours.  Remote  monitoring  simplifies  the 
task  of  measuring  ambient  concentrations 
and  emissions.  Measurements  of  CO  with 
the  OPO  is  possible  at  the  overtone  band 
at  4260  cm* 1  27 .  Many  lines  without 
water  vapor  interference  can  be  found. 
The  absorption  cross  section  is  approxi¬ 


mately  0.4  (atm  cm]*'  **  which  leads  to 
a  measurement  sensitivity  of  0.3  ppmkm. 
The  measurement  of  atmospheric  back¬ 
ground  concentrations  of  CO  in  unpol¬ 
luted  atmospheres  difficult*  but  the 
observation  of  pollution  and  measure¬ 
ments  of  source  emissions  are  possible. 

S02 :  Sulfur  dioxide  is  the  major  pollu¬ 
tant  in  sulfurous  smog.  It  has  been 
recognized  as  a  health  hazard*  and  the 
Clean  Air  Act  requires  annual  average 
SO2  concentrations  of  less  than  0.03 
ppm,  although  over  a  3  hour  period  the 
average  may  be  as  high  as  0.50  ppm.  Ue 
have  demonstrated  the  possiblity  to 
measure  sulfur  dioxide  with  the  OPO 
L1DAR  earlier12.  The  Vi+Vj  combination 
band  around  2499  cm*1  has  1 inestrengths 
of  0.4  (atm  cml'1  2®*  which  makes  meas¬ 
urements  of  concentrations  of  0.25  ppmkm 
possible.  Houever*  this  measurement  has 
to  be  made  near  the  end  of  the  tuning 
range  of  the  Litlb03  .  OPO,  where  the 
energy  output  is  low  and  pulse  energy 
fluctuations  are  increased.  Thus  a 
large  number  of  laser  pulses  have  to 
averaged  to  achieve  the  desired  sensi¬ 
tivity.  However,  measurements  of  source 
emissions  and  monitoring  of  S02  polluted 
atmospheres  is  possible. 

03:  Ozone  is  a  major  constituent  in  pho¬ 
tochemical  smog.  There  is  a  natural 
ozone  background  of  0.01  ppm,  but  its 
concentration  reaches  0.3  ppm  in  photo¬ 
chemical  smog25.  Uithin  the  tuning 
range  of  the  OPO,  the  strongest  absorp¬ 
tion  band  is  the  3V3  band  at  3050  cm*’ 
The  strongest  lines  are  about  0.3  (atm 
cmj"'  ,e*  which  leads  .  to  a  sensitivity 
of  0.3  ppmkm.  This  value  is  not  ade¬ 
quate  to  measure  atmospheric  ozone  with 
the  OPO  L1DAR. 

N0X:  Nitrogen  oxides  are  important  in 
photochemical  smog.  Their  primary 
anthropogenic  source  is  combustion.  In 
polluted  urban  regions  the  concentration 
of  NO*  may  be  0.1  ppm,  or  even  higher  in 
exceptional  circumstances.  The  measure¬ 
ment  of  NO  is  not  possible  with  the  OPO 
IIDAR*  since  only  a  weak  overtone  trans- 


ftion  exists  at  4746  cm*1.  But  NO*  has 
relatively  strong  lines  in  the  V1+V3 
band  at  2910  cm*'.  From  spectra  in 
ref. 29.  the  strongest  line  has  a 
strength  of  3.2  (atm  cm]*',  which  allows 
measurements  with  a  0.03  ppmkm  sensitiv¬ 
ity.  Measurements  of  sources  as  well  as 
monitoring  of  polluted  regions  is  possi¬ 
ble.  However.  potential  interference 
from  water  and  hydrocarbon  lines  makes 
careful  line  selection  necessary  to  min¬ 
imize  interference  with  those  species. 

NjO:  Another  nitrogen  compound  in  the 
atmosphere  is  Nitrous  Oxide.  There  is 
an  atmospheric  N2O  concentration  of  0.29 
ppm.  which  is  mainly  caused  by  biologi¬ 
cal  action  in  soil.  It  is  inert  in  the 
troposhpere  but  of  interest  due  to  its 
interaction  with  ozone  in  the 
stratosphere30.  The  absorption  band  at 
2545  cm" 1  is  strong  enough,  to  promise  a 
sensitivity  of  0.05  ppmkm,  and  interfer¬ 
ence  free  lines  can  be  found10.  This 
allows  measurements  of  atmospheric  back¬ 
ground  concentrations  of  nitrous  oxide. 

NHo s  Another  pollutant  of  importance  is 
Ammonia.  Its  strongest  absorption  bands 
within  the  tuning  range  of  the  0P0  are 
around  3400  cm'1  31 .  However,  we  have 
not  found  spectra  of  adequate  accuracy 
to  determine,  if  absorption  measurements 
are  possible  using  lines  located  in  the 
strong  surrounding  water  vapor  absorp¬ 
tion  region. 

tlCHO :  Formaldehyde  also  is  a  component 
of  photochemical  smog.  Its  atmospheric 
concentration  varies  betueen  1  and  50 
ppb25.  The  cross  sections  of  lines 
around  2800  cm"1  from  ref. 29  promise  a 
sensitivity  of  0.01  ppmkm,  which  makes 
long  path  measurements  suitable  to  moni¬ 
tor  formaldehyde  concentrations  in  smog 
conditions. 

HCI :  The  natural  background  of  Hydrogen 
cloride  is  negligibly  small,  but  anthro¬ 
pogenic  sources  can  release  HCI  in  large 
quantities  that  lead  to  concentrations 
of  1  to  5  ppm25.  HCI  has  very  strong 
absorption  lines  around  2900  cm" 1 .  From 


the  spectrum  in  ref. 29  we  determine 
cross  sections  of  10  [atm  cm]'1,  which 
allows  sensitivities  of  0.01  ppmkm. 
Since  several  lines  can  be  selected, 
interference  from  water  vapor  and  hydro¬ 
carbons  can  be  minimized.  Measurements 
of  HCI  are  possible  with  the  0P0  LIDAR. 

Hydrocarbons:  Hydrocarbons  are  important 
in  photochemical  processes.  Only  meth¬ 
ane  has  a  large  background  in  unpolluted 
atmospheres,  all  other  hydrocarbons 
arise  mainly  from  anthropogenic 
sources25.  The  major  sources  are  com¬ 
bustion  and  chemical  processes.  All 
hydrocarbons  have  a  very  strong  V3 
vibration  band,  but  interferences  with 
water  vapor  lines  as  uell  as  other 
hydrocarbons  cause  problems  for  the 
measurements  of  some  of  them.  Methane 
has  a  background  concentration  in  an 
unpolluted  atmosphere  of  1.5  ppm.  Me 
have  discussed  the  measurement  of  C Hi) 
with  the  0P0  LIDAR  in  section  3  and  dem¬ 
onstrated  the  sensitivity  in  such  a 
measurement. 

The  other  hydrocarbon  components  play  a 
more  important  role  in  photochemical 
smog  than  methane.  However,  due  to  the 
difficulties  to  measure  them  independ¬ 
ently,  little  is  knoun  about  reactions 
of  specific  components  in  the  atmos¬ 
phere.  Mith  the  OPO  LIDAR  a  large  num¬ 
ber  of  hydrocarbons  can  be  measured  with 
good  sensitivity.  However,  uater  vapor, 
methane  and  other  components  all  have  to 
be  measured  simultaneously  due  to  inter¬ 
ference  of  lines.  Using  high  resolution 
spectra  obtained  uith  the  EPA  ROSE 
System32,  we  identified  measurement  pos¬ 
sibilities  for  Ethane,  Propane  and  Ethy¬ 
lene. 

Ethane  can  be  measured  at  2976.8  cm'1, 
where  it  has  an  absorption  line  uith  a 
strength  of  13  [atm  cm]"1.  This  line  is 
overlapped  by  a  water  vapor  line,  but 
transmission  on  this  uavelength  is  ade¬ 
quate  for  the  use  of  long  paths.  This 
line  promises  a  sensitivity  of  8  ppbkm, 
if  errors  due  to  the  measurement  of  H2O 
concentration  can  be  neglected. 

Propane  has  a  very  strong  line  at 
2967.6  cm"1,  the  absorption  cross  sec- 


tion  was  determined  as  83  {atm  eml*’. 
However, Hz0,  CHi,  >  and  other  hydrocarbons 
have  absorption  lines  close  to  this  wav¬ 
elength,  and  their  interference  reduces 
the  sensitivity.  But  a  10  ppbkm  sensi¬ 
tivity  still  can  be  expected. 

Ethylene  has  a  relatively  strong  line 
at  2980  cm*',  but  this  line  is  com¬ 
pletely  obscured  by  a  strong  uater  vapor 
line.  Measurement  is  possible  at  2962.5 
cm"1,  but  linestrength  here  is  only  1.7 
(atm  cm]**,  which  does  not  promise  the 
measurement  of  traces  in  the  atmosphere. 

Other  hydrocarbon  components  have  to  be 
evaluated  in  detail. 


He  demonstated  the  capabilities  of  a 
LiNbOj  OPO  as  a  LIDAR  transmitter.  The 
0P0  has  a  tuning  range  of  4500  cm*1  from 
1.4  to  4.0  pm.  The  linewidth  is  less 
than  1  cm'1,  and  can  be  reduced  to  0.1 
cm*1  by  insertion  of  a  tilted  etalon. 
Wavelength  selection  of  the  OPO  is  done 
under  computer  control.  An  extensive 
software  routine  was  developed  to  con¬ 
trol  automatic  measurements. 

Reliability  of  the  system  was  demon¬ 
strated  in  a  series  of  remote  methane 
measurements  at  3.4  pm  over  a  775  m  long 
atmospheric  path,  using  a  building  as 
non  cooperative  backscattering  target. 
The  longest  of  these  measurements  was 
over  a  period  of  18  hours,  during  which 
the  system  operated  automatically  under 
computer  control.  Our  concentration 
measurements  were  compared  with  those  of 
a  nearby  point  monitoring  station.  The 
agreement  was  excellent. 

The  advantages  of  a  continuously  tunable 
infra  red  transmitter  were  further  dem¬ 
onstrated  with  a  measurement  of  metereo- 
logical  importance.  He  measured  simul¬ 


taneously  average  temperature  and 
humidity  over  the  775  m  long  path.  Me 
selected  for  this  measurement  two  nearby 
lines  in  the  1.9  pm  uater  band,  which 
have  a  1456  cm*’  lower  level  energy  dif¬ 
ference  and  so  promise  a  good 
temperature  sensitivity,  and  are  of  sim¬ 
ilar  strength.  Using  those  lines,  we 
measured  temperature  with  a  1°C  relative 
error,  while  humidity  data  have  a  rela¬ 
tive  error  of  Vi. 

In  addition  to  the  described  measure¬ 
ments,  a  large  number  of  other  trace 
constituents  can  be  measured.  He  used 
available  spectroscopic  data  to  identify 
the  sensitivity  for  measurements  of  the 
most  important  trace  constituents  in  the 
troposphere.  He  found,  that  background 
concentrationes  of  some  pollutants  can 
be  measured.  Source  emission  monitoring 
is  feasable  for  most  important,  pollu¬ 
tants  with  the  OPO  LIDAR. 

The  transmitter  source  can  be  improved. 
An.  increase  of  output  energy  to  more 
than  70  mJ  is  possible.  With  this 
amount  of  transmitted  energy,  depthre- 
solved  measurements  using  Mie  scattering 
as  distributed  reflector  are  possible. 
This  will  have  advantages  for  source 
monitoring  applications.  For  the  meas¬ 
urement  of  ambient  background  concentra¬ 
tions,  however,  a  monostatic  system  like 
the  one  described  here  promises  the 
highest  sensitivity. 
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CM*  measurement 

Temperature  measurement 

Transmitter 

0P0/0PA 

0P0 

frequency  range 

2917-2922  cm*1 

5650-5652  cm*’ 

Energy  /  pulse 

7  mJ 

5  mJ 

Pulse  length 

10  nsec 

10  nsec 

Pulse  rep.  rate 

10  Hz 

10  Hz 

Lineuidth  (FWHJ1) 

1.4  cm’ ’ 

0.1  cm*’ 

Receiver 

Light  col  lector 

Newtonian  telescope 

Newtonian  telescope 

Diameter 

400  mm 

400  mm 

focal  length 

1.2  m 

1.2  m 

Detector 

InSb  at  77  K 

InSb  at  77  K 

Detector  D* 

Electronics 

1.3*10’ ’ (cm^Hz/H) 

4*lO’°(cm\/Hz/U) 

Impulse  response  time 

300  nsec 

300  nsec 

Table  ts  Parameters  of  the  OPO-LIDAR 


SPECIES 


ADSORPTION 
CROSS  SECTION 
I  atm  cm]*1 


AT 

WAVELENGTH 
I  cm  I  *  1 


PREDICTED 
SENSITIVITY 
(ppm  km] 


ATMOSPHERIC 

CONCENTRATIONS 


CO* 

0.02 

4969 

CO 

0.3 

4285 

SO* 

0.4 

2499 

3.2 

2913-2928 

N*0 

2.3 

2578 

HCHO 

*10. 

2898 

HCI 

*10. 

2822-2955 

CHi| 

>40. 

3030-2910 

Ethane 

13. 

2976.8 

Ethylene 

1.7 

2962.5 

Propane 

>50. 

2967.6 

H20 

many  differ 

%ent  lines 

Temp. 

0.005 

5650,5651 

using  H20 

2.5 

background  330  ppm 

0.3 

in  pollution  30-50  ppm 

0.25 

in  sulfurous  smog  1  ppm 

0.03 

in  photochem. smog  0.1  ppm 

0.05 

background  0.29  ppm 

0.01 

in  photochem. smog  0.05  ppm 

0.01 

1-5  ppm  local ly 

0.005 

1.5  ppm  background 

0.010 

concentration  of 

0.060 

all  hydrocarbons  except  CIU 

0.010 

<  0.25  ppm  in  3  hour  average 

n 

1°C 

1-2JC 

Table  2:  Predicted  measurement  capabilities  tor  the  0P0  L1DAR. 


Fig.1  Schematic  of  the  1.4  to  4.0  pm  tunable  transmitter. 

Fig. 2  Schematic  of  the  electronics  and  data  processing  based  on 
PDP-11/10  minicomputer. 

Fig. 3  4S0  cm*1  scan  uith  the  0P0  LIDAR  around  3.0  pm.  The  lower  trace 

shows  the  measured  transmission,  using  a  building  775  m  away  as 
backscatter ing  target.  The  upper  trace  shows  a  simulation  of  the 
same  transmission  spectrum  using  the  AFCRL  data  base.  A  lineu- 
idth  of  1.5  cm*1  was  assumed  for  the  simulation.  But  the  rectan¬ 
gular  lineshape  in  the  simulation  does  not  fit  the  OPO  lineshape 
very  well. 

Fig. 4  Photo  acoustic  spectrum  of  2vo  CHi,  Q  branch  at  1.66  pm..  The  scan 
is  repeated  with  and  without  the  tilted  etalon  in  the  OPO  cavity. 
The  linewidth  of  the  OPO  in  each  case  is  apparent. 

Fig. 5  Transmission  spectrum  around  3.4  pm.  Upper  trace  shows  the 
atmospheric  transmission  over  a  1.55  km  long  path.  The  lower 
trace  is  transmission  through  a  1.2  m  long  absorption  cell  filled 
with  CHi)  and  Hj  and  is  used  to  identify  CHi,  lines  in  the  atmos¬ 
pheric  spectrum. 

Fig. 6  One  hour  segment  of  remote  methane  measurement,  shouing  CHi,  con¬ 
centration  vs.  time.  Each  point  represents  a  20  shot  or  4  sec 
average,  and  has  a  rms~f luctuation  of  0.6  ppm.  The  heavy  line 
shows  same  data  after  later  averaging,  reducing  rms-erros  to  0.09 
ppm. 

Fig. 7  CH*  concentration  vs.  time  over  a  17  hour  period.  Each  point  is 
an  average  over  80  sec,  reducing  rms-deviations  to  0.09  ppm, 
which  is  less  than  the  measured  CHh  fluctuation  in  the'atmos- 
phere.  Open  circles  indicate  measured  concentrations  in  Reduood 
City  point  monitoring  station. 

/* 

Fig. 8  CHi,  concentration  vs.  time.  Stanford  LIDAR  data  are  compared  to 
values  from  Redwood  City  point  monitoring  station  10  km  away. 
Here  a  prevailing  wind  yielded  a  30  min  delay  between  the  atmos¬ 
phere  at  point  monitoring  station  and  its  arrival  in  Stanford. 

% 

Fig. 9  High  resolution  atmospheric  scan  over  1.55  km  path.'  The  OPO 
linewidth  is  0.1  cm'1.  The  OPO  transmission  scan  (lower  trace) 
is  compared  to  a  simulation  using  AFCRL  tapes  (upper  trace). 

Fig. 10  Temperature  and  absolute  humidity  vs.  time.  Measurement  from 
morning  of  11th  march  1980.  Temperature  data  are  compared  uith 
record  from  thermograph  at  receiving 'telescope. 


ABSORPTION  CELL 


OPO  0PA 


RELATIVE  INTENSITY  TRANSMISSION 


ATMOSPHERIC  TRANSMISSION 
OVER  1.55  km  LONG  PATH 


2900  3000  3100  3200 


FREQUENCY  (cm'1) 

Fig.  3*  *»50  cm  scan  with  the  OPO  HOAR  around  3.0  pm.  The  lower  trace  shows  the 

measured  transmission,  using  a  building  775  m  away  as  backseat te r i ng  target. 
The  upper  trace  shows  a  simulation  of  the  same  transmission  spectrum  using 
the  AFCRL  data  base.  A  linewidth  of  1.5  cm"1  was  assumed  for  the  simulation. 
But  the  rectangular  iineshape  in  the  simulation  does  not  fit  the  OPO  lineshape 
very  well. 


OPO  RESOLUTION 
MEASURED  ON 
CH4  Q  BRANCH  IN 
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Photoacoustic  spectrum  of  2vj  CH4  Q  branch  at  1.66  ym.  The  scan  is  repeated  with 
and  without  the  tilted  etalon  in  the  OPO  cavity.  The  llnewldth  of  the  OPO  in  each 
case  is  apparent. 
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Fig.  5.  Transmission  spectrum  around  J.k  ym.  Upper  trace  shows  the  atmospheric 
transmission  over  a  1.55  km  long  path.  The  tower  trace  is  transmission 
through  a  1.2  m  long  absorption  cel)  filled  with  CH4  and  N2  and  is  used  to 
identify  CH4  lines  in  the  atmospheric  spectrum. 
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Mg. 6  One  hour  segment  of  remoto  methane  measurement*  showing  CII*  con¬ 
centration  vs.  time.  Each  point  represents  a  20  shot  or  4  sec 
average,  and  has  a  rms-f luctuation  of  0.6  ppm.  The  heavy  line 
shows  same  data  after  later  averaging,  reducing  rms-erros  to  0.09 


oicof'-t0in<j-fooj--o 


(Ujdd)  N0I1VUXN3DN0D 


Fig. 7  Clt*  concentration  vs.  time  over  a  17  hour  period.  Each  point  is 
an  average  over  80  sec*  reducing  rms-deviations  to  0.09  pp^, 
which  is  Jess  than  the  measured  CH»  fluctuation  in  the  almos- 
phere.  Open  circles  indicate  measured  concentrations  in  Rcduood 
City  point  monitoring  station. 
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flj.4  Clin  concentration  vs.  time.  Stanford  LIDAR  data  are  compared  to 
values  from  Redwood  City  point  monitoring  station  10  km  away. 
Here  a  prevailing  wind  yielded  a  30  min  delay  between  the  atmos¬ 
phere  at  point  monitoring  station  and  its  arrival  in  Stanford. 
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ATMOSPHERIC  TRANSMISSION  OVER  1.55km  LONG  PATH 


High  resolution  atmospheric  scan  over  1.55  km  path.  The  0P0  linewidth  is  0.1  ci 
The  OP 0  transmission  scan  (lower  trace)  is  compared  to  a  simulation  using  AFCRL 
tapes  (upper  trace). 
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rig. 10  Temperature  and  absolute  humidity  vs.  time.  Measurement  from 
morning  of  11th  march  1900.  Temperature  data  arc  compared  tilth 
record  from  thermograph  at  receiving  telescope. 
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APPENDIX  B 


GROWTH  OP  HIGH  PHASEMATCHED  TEMPERATURE  LiNb03  SINGLE  CRYSTALS 
AND  APPLICATION  TO  EFFICIENT  SECOND  HARMONIC 
GENERATION  OF  Nd : YAG 

Robert  L.  Byer  and  Yong  Kwan  Park 
Edward  L.  Ginzton  Laboratory 

and 

Robert  S.  Feigelson  and  Wayne  L.  Kway 
Center  for  Materials  Research 
Stanford  University 
Stanford,  California  94305 


I.  INTRODUCTION 

Optically  induced  refractive  index  damage  in  LiNbC>3  has  been 
known  for. over  ten  years.1  When  LiNbO^  crystals  are  used  in  non¬ 
linear  devices  involving  second  harmonic  generation  (SHG)  refrac¬ 
tive  index  changes  are  induced  which  greatly  reduce  the  conversion 
efficiency  and,  therefore,  the  usefulness  of  the  material.  The 
large  nonlinearity  of  LiNb03  along  with  its  ability  to  achieve 

non-critical  phasematched  SHG  of  the  Nd:YAG  laser  wavelength  at 
2 

1.064  ym  has  led  to  considerable  effort  to  overcome  the  induced 
refractive  index  damage. 

It  was  experimentally  discovered  that  the  optical  index  damage 
in  LiNb03  crystals  is  annealed  at  temperatures  above  approximately 
170  -  200°C. lf ^  5  Therefore,  if  LiNb03  is  grown  such  that  the 
phasematching  temperature  is  above  the  annealing  temperature  the 
index  of  refraction  damage  is  essentially  eliminated  and  efficient 
SHG  of  Nd:YAG  is  possible. 


6  7 

Fay  et.al.,  and  Bergman  et.al.,  reported  that  the  phase¬ 
matching  temperagure  of  LiNbOg  could  be  Increased  by  Increasing 
the  LigO  concentration  in  the  melt.  They  showed  that  the  SHG 
phasematching  temperature  for  1.06  urn  could  be  increased  from 
4°C  for  congruently  grown  crystals  with  48.6  mole  %  Lig 6  to  greater 
than  180°C  for  54  mole  %  LigO.  However,  it  was  also  found  that 
the  Li gp  concentration  and  thus  SHG  phasematching  temperature  in 
the  crystal  varied  along  the  crystal  growth  axis  for  growth 

e 

conditions  away  from  the  congruent  melt  composition.  Congruently 

8 

grown  crystals  do  not  exhibit  this  variation. 

The  cause  of  optical  index  damage  in  LiNbOg  has  been  extensively 

studied  not  only  because  it  is  detrimental  for  nonlinear  optical 
9-13 

applications,  but  also  because  it  can  be  used  as  a  basis 

for  volume  holographic  storage;4,  14  '  ^ 

4 

Chen  porposed  that  the  index  damage  was  due  to  a  modulation 
of  the  index  of  refraction  via  the  electro-optic  effect.  The 
electric  field  was  generated  by  an  internal  space-charge  due  to 
the  drift  of  photo-excited  electrons  from  the  region  of  high  light 
Intensity  to  trap  sites  in  the  area  of  low  light  intensity.  Johnston*7 
later  modified  Chen's  model  by  proposing  that  the  internal  field 
In  the  crystal  could  be  due  to  the  change  in  spontaneous  polar¬ 
ization  in  the  illuminated  region  as  a  result  of  charged  lattice 


defects  and  that  the  divergence  of  the  polarization  was  the  source 
of  the  electric  field.  The  defect  in  LiNbOg  was  identified  as  Fe^+ 
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inpurity  ions  present  in  nominally  pure  crystals  and  the 

2+  3+ 

optical  index  damage  was  attributed  to  the  Fe  to  Fe 
reduct ion . 

This  model  led  to  efforts  to  eliminate  iron  impurities 

in  LiNbOg  and  thus  to  eliminate  or  reduce  the  magnitude  of  the 

optical  index  damage.  It  was  found,  however,  that  even  very  pure 

LINbOg  crystals  showed  optical  damage.  In  an  extension  of  the 
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above  model,  Zverev  et.al.,  *  suggested  that  the  reduction  of 
niobium  in  the  crystal  could  account  for  residual  index  of 
refraction  damage  for  even  pure  crystals.  However,  it  was 
expected  that  the  annealing  temperature  should  be-  lower  due  to 
the  smaller  trapping  energy  of  the  reduced  niobium  site  compared 
to  the  iron  site. 

The  efforts  to  reduce  or  eliminate  the  index  damage  effect 

while  maintaining  crystal  optical  quality  included  both  the 

•  18 

reduction  of  the  iron  impurity  concentration  and  the  growth 

•  a  • 

of  the  high  phasematching  temperature  LiNbOg  (hot  phasematching 

23 

LiNbO.,)  from  special  melt,  compositions,  Bridenbaugh  et.al., 

0  *  “  . 

showed  that  MgO  could  be  used  as  a  co-dopant  along  with  LigO  to 

raise  the  phasematching  temperature  and  simultaneously  to  compensate 

for  excess  LigO  segregation.  However,  it  was  found  that  the  effect 

of  MgO  appeared  to  saturate  near  1.0  mole  %.  Another  suggested 

approach  was  to  grow  LiNbOg  at  the  congruent  melting  point  in  order 

to  maintain  crystal  homogeniety  and  then  to  diffuse  lithium  into 

24-26 

the  single  crystal  to  increase  the  phasematching  temperature. 
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After  a  careful  review  of  the  earlier  LiNbOg  crystal 
growth  studies,  we  initiated  a  program  to  grow  hot  phasematched 
LiNbOg  crystals  with  high  optical  quality  for  efficient  SHG  of 
the  1.064  pm  Nd:YAG  wavelength.  The  program  had  the  following 
objectives:  i)  to  grow  LiNbOg  with  excess  LigO  co-doped  with  MgO 
to  arrive  at  a  composition  that  yielded  high  optical  quality 
crystals  with  a  phasematching  temperature  in  excess  of  180°C  ; 
ii)  to  vary  the  growth  parameters  including  temperature  gradients, 
pull  rate,  and  rotation  rate  to  obtain  optimum  crystal  quality; 
ill)  to  re-examine  the  optical  damage  annealing  temperature  to  re¬ 
define  the  minimum  allowable  phasematching  temperature;  iv)  to  test 
crystals  to  correlate  SHG  nonlinear  conversion  efficiency  with 
growth  parameters;  and  v)  to  optimize  and  define  an  optimum  growth 
process  for  the  production  of  hot  phasematching  LiNbOg  crystals 
for  device  applications. 

II.  CRYSTAL  GROWTH  PROCEDURES 

Crystal  growth  was  carried  out  by  the  standard  Czochralski 
technique  for  LiNbOg  utilizing  an  RF  induction  furnace.  The 
furnace  set  up  is  shown  in  Fig.  1.  The  platinum  crucible"  used 
was  either  2"  dia.  by  2"  high  or  2£"  in  dia.  by  2J"  in  height. 

The  charge  size  ranged  between  259  g  and  280  g  for  the  former  and 
510  g  to  580  g, for.  the  latter.  The  growth  parameters  which 
differed  from  those  fpr  standard  LiNbOg  were  mainly  the  pull  rate 
and  the  rotation  rate. 


-  48  - 


For  the  present  testing  and  immediate  device  application 
the  crystal  samples  were  fabricated  with  cross-sections  of  5  mm 
square  and  with  lengths  of  1  cm  minimum  oriented  with  the  long 
axis  perpendicular  .to  the  crystal  optic  axis.  This  means  that 
the  crystals  could  be  grown  in  either  the  a-,  or  the  c-  direction. 
However,  the  a-axis  grown  crystal  proved  to  be  less  desirable 
because  of  larger  compositional  variations  and  resultant  lower 
optical  quality.  Thus  the  c-axis  was  chosen  as  the  standard 
growth  axis  in  this  investigation. 

The  selection  of  temperature  gradients,  pull  rate,  and 
crystal  rotation  rates  is  commonly  governed  by  the  melt  composition 
and  growth  interface  morphology.  Since  the  crystals  were  grown 
from  melts  containing  excess  LigO,  interface  partitioning  and 
build-up  was  expected.  To  prevent  breakdown  and  constitutional 
supercooling,  both  axial  and  radial  temperature  were  interested 
and  the  growth  rates  used  were  decreased.  For  the  present  system, 
the  temperature  gradients  necessary  for  stable  growth  were  regulated 
by  both  the  RF  coil  position  and  an  after-heater  where  the 
temperature  was  set  at  800°C,  such  that  the  temperature  increases 

below  the  melt  surface  and  decreases  along  the  radial  direction  in 

*•  1 

the  melt  surface.  The  pull  rates  used  were  from  1.6 -4.3  mm/hr,  which 
resulted  in  growth  rates  of  approximately  1.8  -5.0  mm/hr.  The  crystal 
rotation  rates  were  varied  between  10  and  50  rpm.  Since  varying 
the  pull  rate  and  rotation  rate  modified  the  interface  shape  and 
thus  the  temperature  gradient  and  compositional  characteristics 


of  the  crystal  and,  therefore,  the  resultant  optical  quality, 
the  rotation  and  pull  rate  changes  were  made  together. 

Crystal  quality  evaluation  was  based  on  SHG  conversion 
efficiency  tests,  phasematching  temperature,  and  optical  homogeniety. 
These  optical  tests  were  correlated  with  the  interface  shape  which; 
in  turn  was  correlated  with  the  pull  rate  and  crystal  rotation  rate. 

The  optimum  results  were  obtained  on  crystals  with  a  planar 
interface  grown  with  a  pull  rate  of  1.6  mm/hr  and  a  crystal 
rotation  rate  of  50  rpm.  Under  these  conditions  the  SHG  vs 
temperature  phasematching  curve  showed  a  clean  single  peak  with 
a  width  of  0.5°C,  and  a  conversion  efficiency  of  35%.  The  crystal 
contained  essentially  no  striations  or  banding  and  had  an  un- 
dlstorted  isogyre  pattern  indicating  few  optical  inhomogenienties. 

•The  most  important  aspect  in  this  investigation  was  determining 

the  most  appropriate  melt-composition  to  yield  a  high  phasematched 

temperature,  good  optical  quality  and  compositionally  uniform  crystal. 

Raising  the  phasematching  temperature  was  based  on  the  incorporation 

of  48.6  mole  %  LigO.  Crystals  were,  therefore,  grown  from  melts 

containing  LigO  much  in  excess  of  the  congruent  composition. 

Furthermore,  MgO  was  added  to  the  melt  so  that  an  excess  amount 
•  * 
of  L^O  would  be  incorporated  into  the  crystal.* 
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As  shown  by  the  phase  diagram  in  Fig.  2  and  previous  work, 

Li^O  in  excess  of  the  congruent  composition  is  incorporated 

in  the  crystal  according  to  the  equilibrium  distribution  coefficient 

kc(Exc.  Li20)  -  --2— ^  »  0.3  (1) 

^L(Exc .  Li20) 

Of  course,  the  actual  incorporation  is  a  function  of  the  growth 
parameters  involved.  The  segregation  of  excess  LigO  and  melt 
composition  variation  during  growth  and  consequently  the  phase- 
matching'  temperature  variation  along  the  length  of  the  crystal 

is  large  as  governed  by  the  equation  for  normal  freezing, 

♦  * 

•  * 

Cs  *  kCyi-g)1^1  (2) 

where  Cg  is  the  concentration  of  a  component  in  the  solid,  Cq 
the  concentration  in  the  melt,  k  the  distribution  coefficient, 
and  g  the  fraction  solidified.  As  mentioned  before,  the 
addition  of  MgO  to  the  melt,  results  in  an  excess  amount  of  Li20 

over  the  congruent  composition  in  the  drystal.  As  pointed  out 

% 
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by  Bridcnhaugh,  an  addition  of  1  mole  %  MgO  to  the  melt  results 
in  the  initial  LiNbOg  composition  containing  0.6  mole  %  excess  Li20 
Since  the  distribution  coefficient  for  MgO  in  LiNbO^  is  1.3,  it 
is  being  rapidly  depleted  from  the  melt  as  the  crystal  is  grown. 
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Since  the  excess  L120  is  determined  by  the  M^O  concentration, 

i  the  contribution  for  M  O  decreases  in  its  concentration  in  the 

© 

melt  decreases:  The  Li20  concentration  in  the  crystal,  which 

is  due  to  excess  LigO  added  to  the  melt,  increases  as  growth 

proceeds  and  because  of  these  two  compensating  effects,  a  flat 

Li20  profile  along  the  crystal  length  results. 

Hie  amount  of  excess  Li20  and  MgO  to  be  added  to  the  melt 

for  the  desired  phasematching  temperature  hs  based  on  available 

7 

literature  data.  From  the  results  of  Bergman  et.al.,  and 
6 

Fay  et.al.,  the  slope  of  the  phasematching  temperature  vs 

excess  Li20  was  determined.  Together  with  MgO  doping  experiment, 

23  28 

done  at  Stanford  and  those  reported  by  Bridenbaugh,  ’  extra¬ 
polated  curves  were  plotted  to  estimate  the  initial  melt  compositions. 
Crystals  were  grown  from  LigO-rich  melts  at  different  Li20 
concentrations  with  MgO  doping  ranging  from  1.0  to  2.5  mole  %. 

Early  results  showed  that  MgO  did  not  contribute  to  increasing 
the  phasematching  temperature  beyond  1.0  mole  %.  Consequently 
the  MgO  doping  was  kept  at  1.0  mole  %  for  all  experimental  runs 
with  different  LigO  concentrations  .  The  results  were  generally 
predictable  with  respect  to  the  phasematching  temperature  vs 
concentration  and  fitted  closely  the  initially  extrapolated 
curves.  The  proper  Li20  concentration  to  achieve  a  phasematching 
temperature  of  180°C  along  with  a  fixed  1.0  mode  %  MgO  concentration 
was  determined.  The  Li20  concentration  vs  phasematching  temperature 
is  shown  in  Fig.  3. 
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Aside  from  the  growth  problems,  poling  is  difficult  in 
hot  phasematch  crystals.  As  the  phasematching  temperature 
reaches  180°C  or  above,  the  margin  of  error  for  poling  temperature 
is  only  20°C.  At  too  low  a  temperature  the  crystal  will  not  be 
completely  poled  and  if  too  high  the  crystal  will  melt.  For 
these  reasons  it  is  desirable  to  go  to’  a  lower  phasematch  temperature 
where  both  growth  and  poling  can  be  carried  out  with  greater  ease. 

Although  no  attempt  has  been  made  in  the  reduction  or 

+2 

elimination  of  the  Fe  impurity,  the  re-evaluation  of  the  self¬ 
annealing  temperature  was  carried  out  because  it  was  believed  that 
the  LiNbOg  crystals  grown  at  Stanford  have  been  handled  with  great 
care  to  achieve  low  impurity  content.  The  measurements  showed 
that  the  self-annealing  temperature  was  indeed  at  an  unexpected 
low  value  of  about  120°C.  With  this  information,  crystals  were 
grown  with  the  new  phasematching  temperature  target  of  120°C  -  130°C. 

The  crystals  with  phasematching  at  about  120°C  showed  all 
parameters  to  be  as  good  as  or  better  than  the  high  phasematch 
crystals  produced  earlier.  The  growth  procedures  were  thus  finalized 
for  producing  crystals  with  the  target  phasematching  temperature 

of  130°C.  The  following  is  the  summary  of  procedures  for  the 

* 

routine  production  and  testing  of  high  phasematch  LiNbOg  crystals 
for  high  power  laser  applications.  . 

The  starting  materials  were  high  purity  LigCOg 
Nb^Og,  and  MgO.  They  wore  weighed  out  in  the  proper  proportion 
to  yield  the  composition  of  .505  Li20,„.495  NbgOg  and  1.0  mole  %  MgO, 
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for  a  phaseraatching  temperature  of  130°C,  well  mixed, and  melted 
t  down  in  a  platinum  curcible  utilizing  a  low  temperature  gradient 
furnace.  Generally,  the  reaction  and  melt  down  is  carried  out 
in  a  box  type  of  SIC  furnace  heated  at  a  rate  of  approximately 

300°C/hr  to  1 ,280°C  and  soaked  at  that  temperature  for  thirty 

■ 

minutes.  At  the  end  of  the  cycle  the  furnace  was  shut  down  and 
the  material  was  allowed  to  cool  in  the  furnace  to  room  temperature. 
This  procedure  was  repeated  until  the  crucible  was  filled  at  260- 
280  for  the  2"  crucible  and  560  -  580  for  the  2}”  crucible. 

Crystal  growth  was  accomplished  in  an  RF  induction  heated  furnace 

•  \ 

with  the  platinum  curcible  as  susceptor.  The  furnace  was  set  up  ; 

1 

with  the  proper  size  coil  and  alumina  parts  for  either  the  2"  dia. 
by  2"  height  or  2i"  dia.  by  2i"  height  crucible.  Together  with 
an  after-heater  maintained  at  800°C,  the  coil,  crucible,  and 
furnace  parts  were  positioned  and  adjusted  to  produce  positive 

*  t 

temperature  gradients  in  the  furnace  with  radial  gradient  of  +5°C 
from  center  to  edge  at  the  melt  surface.  -Crystal  growth  was 
commenced  with  a  c-axis  seed,  allowed  to  come  out  to  a  proper 
diameter  and  proceeded  with  that  diameter  until  a  desired  length 
was  reached.  During  the  growth  cycle,  the  temperature  was 
controlled  at  a  programmed  cooling  rate  of  5-7°/hr.  for  shouldering 
and  a  rate  of  ,2-l°/hr.  for  the  body  of  the  crystal.  At  the  end  ] 

of  the  growth  cycle,  the  crystal  was  tapered  off  by  heating  rapidly 
at  10  -  200 /hr, pulled  into  the  after-heater,  and  heated  to  1,100°C.  j 


-  54  - 


Vitb  the  main  furnace  power  (RF)  turned  off  slowly  (taking  about 


i  thirty  minutes),  the  crystal  was  cooled  in  the  covered  after¬ 
heater  to  room,  temperature  at  a  rate  of  100°/hr. 

The  poling  was  accomplished  with  a  flat  temperature  profile 
tubular  SiC  furnace.  The  crystal  was  prepared  for  poling  by 
cutting  off  the  ends  and  attaching  the  platinum  electrode  with 
platinum  paste.  The  crystal  was  placed  in  the  furnace  with  leads 
attached,  was  heated  up  to  1,170°C  at  200°/hr,  and  at  that  point, 
the  poling  field  was  applied. The  poling  voltage  was  adjusted  for 
2  V/cra  and  the  temperature  was  maintained  at  1,200°C  for  a  period 
of  two  hours.  At  the  end  of  this  period  the  temperature  was 
lowered  at  a  moderate  rate  from  1,200°C  to  1,100°C  in  fifteen 
minutes.  At  that  point  the  poling  field  was  turned  off  and  the 
furnace  temperature  was  programmed  to  lower  at  a  rate  of  100°/hr 
all  the  way  to  room  temperature.  During  the  entire  poling  cycle 
oxygen  flow  was  maintained.  The  flow  rate  was  3  cfh  from  room 
temperature  to  1,170°C,  2  cfh  for  poling  and  3  cfh  from  1,200°C 
to  room  temperature.  When  the  crystal  was  removed  from  the 
furnace  after  poling,  the  electrodes  were  stripped  ;and  the  ends 
were  re-cut  and  polished.  The  Domain  structure  was  checked  by 

t 

etching  the  crystal  in  a  3  :  L  HNOg/IIF  solution  at  approximately 
70°C  for  thirty  minutes.  By  examing  the  etched  surfaces  under  a 
microscope,  the  ferro-electric  domain  pattern  was.  revealed.  If 
the  crystal  still  contained  multi-domain  structure,  the  poling 
cycle  was  repeated.  In  some  cases  higher  voltage  or  longer  poling 
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time  was  used.  If  the  crystal  was  completely  poled  and  showed 
a  good  distortion-free  optical  isogyre  pattern,  it  was  ready 
for  fabrication  and  testing.  Figure  4  shows  photographs  of 
boules  515  and  528  which  were  grown  and  poled  by  this  procedure. 

III.  CRYSTAL  EVALUATION 

It  was  realized  at  the  outset  of  this  program  that  progress 
in  the  growth  of  hot  phasematching  LiNbOg  crystals  would  depend 
upon  rapid  testing  and  evaluation.  In  addition  to  optical  quality 
evaluation,  each  boule  was  tested  by  SHG  using  a  cw  lamp  pumped, 
aeousto-optically  Q-switched  Nd : YAG  laser.  The  temperature  phase¬ 
matching  curve  was  obtained  and  the  SHG  efficiency  was  measured 
when  appropriate.. 

•  Initially  crystals  were  prepared  with  phasematching  temperatures 

*  •  * 

at  or  above  the  expected  annealing  temperature  range  of  180-200°C. 

Table  1  summarizes  crystal  parameters  for  a  set  of  hot  phasematching 

crystals  initially  tested.  Among  the  crystals  in  Table  I,  crystal 

#515  shows  the  expected  high  phasematching  temperature  and  single 

narrow  phasematching  peak  as  shown  in  Fig.  5a.  The  reasonably  good 

optical  quality  of  crystal  #515  was  also  verified  by  SHG  conversion 

efficiency  results.  However,  the  phasematching  temperature  of 

the  crystal  increased  by  10°C  from  the  seed  to  the  bottom  end  of 

6 

the  boule.This  result  was  consistent  with  an-  earlier  report 
and  implied  an  inhomogeneous  lithium  concentration  along  the 

5  g  7 

crystal  growth  axis.  '  ' 


LiNbO,  Crystal  Parameters  for  the  Phasematching  Temperature  Greater  than  160°C 


Crystal  #528  was  an  attempt  to  duplicate  crystal  #515. 

The  phasematching  temperature  tuning  curve  shown  In  Fig.  5b 
and  the  SHG  conversion  efficiency  of  crystal  #528  showed 
comparable  optical  quality  with  that  of  crystal  #515.  However, 
the  discrepancy  In  phasematching  temperature  from  crystal  #515 
indicated  that  the  duplication  of  hot  phasematching  LiNb03 
crystals  has  not  been  consistently  obtained  mainly  due  to  the 
difficulty  involved  in  growing  these  crystals  from  the  highly 
non-congruent  melt . 

Another  suggested  hot  phasematching  crystal  preparation 
approach  was  the  growth  and  indiffusion  of  LiNbO^.  This  approach 
was  briefly  attempted  for  some  crystals.  The  crystals  were  grown 
at  the  congruent  melting  point  and  then ‘were  diffused  with  55 
mole  %  Li20  mixture  powder  at  1100°C  for  five  or  ten  hours. 

The  crystals  showed  significant  increase  of  the  phasematching 
temperature  but  severe  deterioration  of  the  optical  quality. 

Following  the  initial  growth  of  high- phasematching 
temperature  crystals  a  re-evaluation  of  the  annealing  temperature 
was  undertaken  to  attempt  to  define  an  optimum  phasematching 
condition.  The  annealing  tests  were  accomplished  with  a  high 
average  power  cw  argon  ion  laser  operating  at  0.5145  ym  which  is 
near  the  0.532  ym  second  harmonic  wavelength* of  Nd:YAG 

The  optical  damage  was  observed  in  two  ways;  by  the  distorted 
beam  shape  of  the  transmitted  beam  for  a  preliminary  observation; 
and  by  the  dynamic  light  scattering  pattern  of  the  transmitted 
beam  for  a  quantitative  measurement.  Both  effects  were  caused 
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by  the  optical  index  damage  of  the  crystal  due  to  the  incident 

.5145  pm  beam.  In  the  preliminary  observation  the  distorted 

transmitted  beam  was  imaged  on  Polaroid*  film  at  various  crystal 

4  7 

temperatures.  Consistent  with  the  earlier  reports,  '  we 
have  observed  that  index  damage  occurred  much  more  severely 
on  the  extraordinary  refractive  index  (n  )  than  on  the  ordinary 
refractive  index  (n^).  Thus  the  measurement  of  the  beam  distortion 
due  to  the  index  damage  was  made  mainly  for  n  wave.  As  an 
example.  Fig.  6  shows  the  evolution  of  the  beam  distortion  for 
crystal  #528.  It  demonstrates  that  the  beam  distortion,  in  turn 
optical  index  damage,  is  eliminated  at  the  unexpectedly  low 
crystal  temperature  of  approximately  120°C.  Other  recently  grown 
LiNbO^  crystals  showed  similar  results.  However,  crystals  grown 
pproximately  ten  years  earlier  showed  annealing  temperatures  near 
the  expected  170-180°  range. 

To  measure  the  annealing  temperature  of  the  crystals  more 
accurately  under  high  power  density  conditions,  the  cw  .5145  ym 
beam  was  focused  very  tightly  into  the  crystal  and  the  scattered 
light  pattern  of  the  transmitted  beam  was  observed  on  the  screen. 
The  transmitted  beam  initially  showed  a  dynamic  scattered  light 
pattern  at  low  crystal  temperatures.  An  increase  of  the  crystal 
temperature  changed  the  dynamic  light  scattering  pattern  to  a 
static  light  scattering  pattern  because  of  the  self-annealing 
.of  the  optical  index  damage.  The  temperature  at  which  the  dynamic 
pattern  changed  into-  the  static  pattern  was  defined  to  be  the 
optical  index  damage  annealing  temperature.  The  annealing 
temperatures  were  measured  for  both  the  0.488  and  0.5145  ym 
argon  ion  laser  lines.  The  annealing  temperature 
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difference  was  negligibly  small  as  observed  in 

previous  reports.  Since  the  index  damage  is  reduced  for  longer 
«  29 

wavelengths  it  is  reasonable  to  consider  that  the  annealing 
temperature  at  .5145  pm  and  at  .5320  pm  are  similar. 

Figure  7  shows  a  plot  of  the  annealing  temperature  vs 
incident  power  density  for  a  hot  phasematching  LiNbOg  crystal 
#528  and  a  congruently  grown  LiNbOg  crystal.  Other  hot  phasematching 
LiNbOg  crystals  showed  the  similar  results.  Although  the  required 
annealing  temperature  of  the  hot  LiNbOg  crystal  is  increased  as 

2 

the  incident  power  density  increases,  it  saturates  at  a  0.1  MW/cm 
power  density  and  eventually  clamps  strongly  below  approximately 

O 

110°C  up  to  1.2  MW/cm  power  density.  Furthermore,  the  extra¬ 
polation  of  the  plot  shows  that  the  crystal  temperature  of  120°C 
should  be  high  enough  to  anneal  the  index  damage  even  at  higher 
power  densities.  Therefore,  contrary  to  the  earlier  LiNbOg  crystals 
which  showed  annealing  temperatures  near  170°C,  the  recently 
Stanford  grown  hot  phasematching  LiNbOg  crystals  show  lower 

annealing  temperatures  of  120°C.  The  low  annealing  temperature 

2+ 

is  consistent  with  an  extremely  low  Fe  impurity  concentration  of 

16—20 

1  ~  2  ppm  in  the  Stanford  crystals. 

The  discovery  of  the  low  120°C  annealing  temperature 

* 

immediately  led  to  the  growth  of  crystals  with  phasematching 
* 

temperatures  just  above  that  value.  The  lower  annealing  temperature 
allowed  a  reduction  of  the  excess  lithium  in  the  melt  and  thereby 
greatly  lessened  the  growth  difficulty  and  enhanced,  the  possibility 
for  growth  of  larger  high  optical  quality  crystals  needed  for 
efficient  SHG. 
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Table  II  lists  the  good  optical  quality  hot  phasematching 


LiNbO,  crystals  grown  at  Stanford  whose  phasematching  temperatures 

w 

were  lower  but  still  high  enough  to  anneal  the  optical  index 

damage  produced  by  the  generated  second  harmonic  of  the  1.06  ym 

beam  at  approximately  20  MW/cm  peak  power  density.  The  second 

harmonic  phasematching  curves  of  crystals  #557  and  #558  are  shown  - 

In  Fig.  8  .  Among  the  crystals  listed,  crystal  #557  shows  the 

lowest  phasematching  temperature.  It  also  shows  that  the  inhomogeneity 

of  the  lithium  concentration  along  the  crystal  growth  axis  is 

greatly  reduced  since  the  phasematching  temperature  variance 

from  the  seed  to  the  bottom  end  of  the  boule  is  decreased  from 

10°C  for  the  previously  grown  hot  LiNbO^  crystals  to  2°C  for 

crystal  #557.  Crystal  #558  is  an  attempt  to  duplicate  the 

crystal  #557  and  shows  that  a  consistent  duplication  with  high 

optical  quality  is  indeed  possible  at  this  melt  composition. 

The  2  -  3°  C  phasematching  temperature  difference  between  two 
crystals  is  mainly  due  to  the  SHG  temperature  measurement  error 
of  the  crystal  oven  controller. 

In  addition  to  complete  index  damage  annealing  and  the  low 
L^O  composition  inhomogeneity,  the  high  optical  quality  of 
these  crystals  is  supported  by  high  SHG  conversion  efficiency. 

The  SHG  conversion  efficiency  of  the  crystals  was  measured  by 
both  TEMco  mode  and  multi-transverse  mode  beam  from  the  1.06  ym 
Nd.'YAG  source.  When  a  TEMqo  mode  beam  with  3  watts  average  power 
and  1  Kllz  pulsed  Q-switch  repetition  rate  was  focussed  into  the 
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2 

2  cm  length  crystals  at  approximately  20  MW/cm  peak  power 
density,  the  crystals  showed  a  40%  average  power  conversion 
efficiency. 

A  10  mrad  divergence  multi-transverse  mode  beam  with  the 

same  peak  power  density  level  produced  a  33%  average  power 

conversion  efficiency.  The  acceptance  angle  of  a  2  cm  long 

31 

LiNbOg  crystal  is  calculated  to  be  15  mrad.  The  comparable 

SHG  conversion  efficiency  of  the  crystals  for  the  multi-transverse 

mode  to  that  for  the  TEM  mode  verifies  that  advantage  of  non- 

oo 

critical  phasematching  is  fully  utilized. 

The  measurements  of  phasematching  temperature  tuning  curves, 
annealing  temperatures,  phasematching  temperature  variance  along 
the  crystal  growth  axis,  and  SHG  conversion  efficiency  demonstrates 
that  high  quality  hot  phasematching  LiNbOg  crystals  can  be  grown 
with  consistent  SHG  performance. 

IV.  CONCLUSION 

During  the  course  of  this  program,  we  have  demonstrated  that 
hot  phasematching  LiNbOg  crystals  can  be  grown  with  the  optical 
quality,  reproducibility  and  size  required  for  device  applications 
by  the  proper  choice  of  composition  and  growth  parameters. 

More  significantly  and  specif icially,  we  have  demonstrated 
that  LiNbOg  crystals  with  a  phasematching  temperature  greater 
than  180°C  can  be  grown  with  small  concentration  gradients  from 
a  melt  doped  with  MgO  and  excess  LigO  together  with  a  planar 
growth  interface;  that  the  phasematch  temperature  can  be  readily 
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adjusted  simply  by  changing  the  excess  LigO  concentration 

i  proportional  to  the  phasematch  temperature  desired;  that  the 

crystal  quality  and  SHG  characteristics  can  be  maintained  as 

long  as  the  MgO  concentration  and  growth  parameters  are  held 

constant.  We  have  also  found  that  the  self-annealing  temperature 

is  much  lower,  near  120°C  for  the  Stanford  grown  crystals,  as 

compared  to  the  previously  reported  value  of  170  -  180°C. 

Taking  advantage  of  the  lower  self-annealing  temperature, 

the  crystals  grown  with  phasematching  temperature  of  approximately 

120°C  showed  exceptional  optical  homogeneity,  and  a  40%  average 

power  SHG  conversion  efficiency. 

The  improved  hot  phasematching  LiNbOg  crystals  along  with 

the  electro-optical  or  piezo-optical  compensation  for  self 
32 

heating  effects  should  greatly  enhance  crystal  device 
applications  to  high  average  power  second  harmonic  generation  of 
the  Nd:YAG  source. 
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FIGURE  6APTI0NS 


A  Schematic  of  the  LiNbOg  crystal  growth  furnace. 

1.  Alumina  Pull  Rod,  2.  After  Heater,  3.  TC  Signal 
to  Aux.  Controller,  4.  Pt.Seed  Holder,  5.  LiNbOg  Crystal 

6.  LiNbOg  Melt,  7.  Pt.  Crucible,  8.  Alumina  Plate, 

9.  RF  Coil,  10.  Alumina  tubes,  11.  Firebricks. 

12.  Sapphire  Rod-Radiation  Pyrometer,  13.  Temp.  Signal 

to  Main  Controller,  14.  Fiberfrax. 

Phase  diagram  of  the  Li20-Nb205  System  (after  Carruthers 
et.al.5). 

Phasematching  temperature  vs  concentration  of  Li20  and 
MgO  in  the  melt. 

Photographs  of  high  phasematching  temperature  LiNbOg 
crystals  (#515  and  #528).  The  crystal  #528  is  an  attempt 
to  duplicate  #515.  Note  the  solid  liquid  interface  shape 
at  the  bottom  of  the  boules! 

Phasematching  temperature  tuning  curve  of  the  crystals 
from  boule  #515  and  #528.  Both  were  sampled  at  the  bottom 
of  the  boules. 

Photographs  of  the  0.5145  pm  beam  propagating  as  an 
extraordinary  wave  through  crystal  #528.  The  index  damage 
is  annealed  at  temperatures  greater  than  120°C. 


Annealing  temperature  vs  power  density  for  0.5145  ym 
incident  on  recently  grown  LiNbO^  crystals. 

Phasematching  temperature  tuning  curve  of  the  crystals 
from  boules  #557  and  #558.  Upper  traces  are  sampled  near 
the  seed  of  the  boules  and  lower  ones  from  the  bottom. 
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PHASE  MATCHING  TEMPERATURE  (#C) 


MOLE  FRACTION  OF  LigO  IN  MELT 

FIG.  3— •Phasematching  temperature  vs  concentration  of  LiO,  and  MgO  in  the 
melt. 
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FIG.  4— Photographs  of  high  phasematching  temperature  LiNb03  crystals  (#515  &  #528).  The  cryatal 
#528  Is  an  attempt  to  duplicate  #515.  Mete  the  solid  liquid  Interface  shape  at  the 
bottom  of  the  boules. 


FIC.  6— -Photographs  of  the  0.5145  pm  beam  propagating  as  an  extraordinary 
wave  through  crystal  0 528.  The  index  damage  is  annealed  at 
tc«rf>craturcs  greater  than  120°C. 
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#558 


•Fhasematchlng  temperature  tuning  curve  of  the  crystals  from  boule  1557  &  #558.  Upper  traces  are 
sampled  near  the  seed  of  the  boules  and  lower  ones  from  the  bottom. 


